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ABSTRACT 

We present results from a survey designed to probe the star formation properties of 32 damped 
Lya systems (DLAs) at z ~ 2.7. By using the “double-DLA” technique that eliminates 
the glare of the bright background quasars, we directly measure the rest-frame FUV flux 
from DLAs and their neighbouring galaxies. At the position of the absorbing gas, we place 
stringent constraints on the unobscured star formation rates (SFRs) of DLAs to 2<r limits of 
Ip < 0.09— 0.27 Mq yr“^,con‘esponding to SFR surface densities Egfr < 10“^ ®— 10“^ ® Mq 
yr“^ kpc“^. The implications of these limits for the star formation law, metal enrichment, and 
cooling rates of DLAs are examined. By studying the distribution of impact parameters as a 
function of SFRs for all the galaxies detected around these DLAs, we place new direct con- 
straints on the bright end of the UV luminosity function of DLA hosts. We find that ^ 13% 
of the hosts have tp ^ 2Mq yr“^ at impact parameters 6dia ^ {iP/Mq yr“^)° ® + 6 kpc, 
differently from cun'ent samples of confirmed DLA galaxies. Our observations also disfavor 
a scenario in which the majority of DLAs arise from bright LBGs at distances 20 ^ 6dia < 
100 kpc. These new findings corroborate a picture in which DLAs do not originate from 
highly star forming systems that are coincident with the absorbers, and instead suggest that 
DLAs are associated with faint, possibly isolated, star-forming galaxies. Potential shortcom- 
ings of this scenario and future strategies for further investigation are discussed. 

Key words: Stars: formation - galaxies: evolution - galaxies: high-redshift - quasars: ab- 
sorption lines - ultraviolet: galaxies - ISM: atoms 


1 INTRODUCTION 

The continuous improvement of the instrumentation available at 
ground based observatories and on board of the Hubble Space Tele- 
scope (HST) has facilitated numerous deep galaxy surveys across 
wide areas of the sky to chai'acterize the typical star' formation 
rates (SFRs), dust content, stellar masses, and structural proper- 
ties of galaxies across cosmic time, up to the epoch of reioniza- 
tion. However, the redshift-dependent sensitivity of these surveys 
offers only a view of the tip of the iceberg of the galaxy popula- 
tion in the distant Universe. An appealing alternative to explore the 
physical properties of galaxies across a wider range of luminos- 
ity and masses is the study of damped Lyman-a absorbers (DLAs; 
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Wolfe et al. 1986, 2005) that are detected along the line of sight 
to bright background sources such as quasars or the afterglows of 
7 -ray bursts (GRBs). Being selected solely according to their neu- 
tral hydrogen column density (Ae 2 X 10^° cm ^), DLAs arise 
in dense and neutral gas and, as suggested also by models and hy- 
drodynamic simulations, they appear to trace the hydrogen content 
of galaxies over a wide range of halo masses (e.g. Nagamine et al. 
2007; Pontzen et al. 2008; Tescari et al. 2009; Hong et al. 2010; 
Cen 2012; Bird et al. 2014). 

The study of hydrogen absorption lines in lai'ge spectroscopic 
samples of DLAs (e.g. Prochaska et al. 2005; Prochaska & Wolfe 
2009; Noterdaeme et al. 2009, 2012a; Zafar et al. 2013) and of the 
associated metal lines in high-resolution spectra (e.g Ledoux et al. 
2006; Prochaska et al. 2007; Rafelski et al. 2012; Jorgenson et al. 
2013; Mplleret al. 2013; Neeleman et al. 2013; Rafelski et al. 2014) 
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have constrained the kinematics and metallicities of the gas in these 
absorbers at different redshifts. DLAs at redshifts z = 2 — 3 cover 
roughly one third of the sky and evolve in time such that the shape 
of the column density distribution function is preserved (Zwaan et 
al. 2005; Prochaska & Wolfe 2009; but see Braun 2012). They are 
typically sub-solar in composition, but more enriched (~ 5% of the 
solar metallicity) than the intergalactic medium (IGM; Pettini et al. 
1997; Prochaska et al. 2007; Fumagalli et al. 2011; Rafelski et al. 
2012; Jorgenson et al. 2013). Finally, DLAs exhibit abundance pat- 
terns that can be described with bursty star formation histories that 
are common in irregular and dwarf galaxies in stochastic regimes 
(Dessauges-Zavadsky et al. 2007; da Silva et al. 2014; Dominguez 
et al. 2014), and, as suggested by recent studies, that follow the 
chemical composition observed in Galactic halo stars (Rafelski et 
al. 2012) or local dwarfs (Cooke et al. 2014). 

Flowever, an outstanding challenge that has limited our under- 
standing of the nature of DLAs, and their link to the galaxy pop- 
ulation at high redshift, is the identification of the host galaxies 
of the absorbers. In fact, the bright background light that enables 
absorption spectroscopy becomes a limiting source of contamina- 
tion during imaging observations, hampering the identification of 
galaxies at close impact parameters to the quasars (e.g. Moller & 
Warren 1998). For this reason, despite over twenty years of inves- 
tigation (see Appendix B of Fumagalli et al. 2010b), the emitting 
counteiparts of only a dozen DLAs have been detected to date at 
z > 2 (Table 2; see also Krogager et al. 2012). Many of these de- 
tections have stemmed from recent specb'oscopic and integral field 
unit (IFU) searches (Fynbo et al. 2010, 2011; Peroux et al. 2012; 
Noterdaeme et al. 2012b; Fynbo et al. 2013; Bouche et al. 2013; 
Jorgenson & Wolfe 2014; Rafelski et al. 2014). 

Despite these recent successes, it is difficult to assess the ex- 
tent to which the identified DLA counterparts represent the gen- 
eral population of DLA host galaxies; this is due to both the se- 
lection criteria used in some searches (e.g. targeting high metallic- 
ity DLAs) and the absence of a reliable census of non-detections. 
Furthermore, by design, the bright background source used for ab- 
sorption spectroscopy is perfectly aligned to the absorbing mate- 
rial, making a direct study of the in-situ SFRs of DLAs detected 
against quasars impossible even in surveys that target representa- 
tive samples of DLAs at the high resolution offered by HST (e.g. 
Warren et al. 2001). In turn, this prevents a systematic analysis of 
the star formation law in neutral gas clouds using far ultraviolet 
(FUV) tracers in the distant Universe. And, in fact, previous work 
on the subject has been limited to Lya emission (Rahmani et al. 
2010; Noterdaeme et al. 2014; Cai et al. 2014) or indirect methods 
only, such as statistically connecting DLAs to low surface bright- 
ness galaxies (Wolfe & Chen 2006) and to the outskirts of compact 
Lyman break galaxies (LBGs; Rafelski et al. 201 1). 

Pressing questions on what the typical DLA counterparts are 
or at what rate these H I clouds form stars remain open. Are DLA 
galaxies massive rotating disks as suggested by the absorption line 
profiles (e.g. Prochaska & Wolfe 1997) or the small building blocks 
predicted by the theories of the hierarchical assembly of structures 
(e.g. Haehnelt et al. 1998; Rauch et al. 2008)? Or are DLAs gas 
clouds in the halo of star-forming galaxies (e.g. Bouche et al. 2013), 
possibly hosting vigorous starburst-driven outflows (e.g. Noter- 
daeme et al. 2012b)? And what is the source of heating that bal- 
ances the cooling rates estimated from the [C It*] absorption lines 
(Wolfe et al. 2003b,a)? 

Recently, two promising methods to address some of the 
above questions have been explored. The first approach focuses on 
intervening DLAs detected along the line of sight to GRBs, taking 


advantage of the fact that once the bright afterglow has faded, deep 
imaging follow-up can be conducted to search for DLA counter- 
parts in emission at all impact parameters. The downside of this 
approach is that the transient nature of GRBs limits the size of the 
samples for which good-quality high-resolution spectra are avail- 
able to study the DLA properties in absoiption (e.g. Cucchiara et 
al. 2014). At present, only a handful of GRB sightlines with in- 
tervening DLAs have been searched for the optical counteiparts of 
these absorbers (Schulze et al. 2012). 

The second approach is to exploit the fortuitous alignment of 
two optically-thick absorbers along individual quasars (Steidel & 
Hamilton 1992; O’Meara et al. 2006; Christensen et al. 2009). The 
basic idea of this technique, which we dub the “double-DLA” or the 
“Lyman limit” technique, is to image quasar fields that host DLAs 
at lower redshift and optically-thick absorbers (DLAs or LLSs, that 
is Lyman limit systems) at higher redshift. As discussed in Fuma- 
galli et al. (2010b) and Fumagalli et al. (2014b), this second higher- 
redshift system acts as a natural blocking filter that absorbs all the 
quasar light blueward of its Lyman limit, potentially allowing the 
detection of DLA counterparts at all impact parameters at these 
shorter wavelengths. Nevertheless, given the high intrinsic fluxes of 
the background quasars, positive detections should be carefully vet- 
ted to exclude residual contamination from the background sources 
(see Cai et al. 2014, Section 2.1.1). The main limitation of this ap- 
proach is that it requires a specific alignment of two absorbers both 
in space and redshift. However, thanks to the large volume probed 
by modern spectroscopic surveys (e.g. Noterdaeme et al. 2012a), 
this experiment can be performed for several quasars with inter- 
vening DLAs at z > 2. 

In previous papers of this series (Fumagalli et al. 2010b, 
2014b), we discuss in detail the survey rationale, and we demon- 
strate the power of this double-DLA technique in avoiding the con- 
tamination of the bright background quasars, which enables us to 
achieve the same sensitivity limit at all impact parameters from the 
DLAs. In this previous work, we also describe the imaging and 
spectroscopic observations of the 32 quasar fields with intervening 
DLAs at z ~ 1.9 — 3.8, which comprise our sample. Information 
on the data reduction, survey completeness limits, and preparation 
of the galaxy catalogues can be found in Fumagalli et al. (2014b), 
where we also present the Ah i and metallicity distribution of our 
sample, as measured in absorption with echellette spectra. In our 
previous work, we also show that the targeted DLAs constitute an 
unbiased sample that represents the full range of absorption prop- 
erties of the parent DLA population at these redshifts. 

In this third paper, we discuss the results of this survey. In 
Section 2 we derive the first direct limits for the observed in-situ 
SFRs of DLAs in the UV, also discussing the implications for the 
metal content and the heating rates of DLA gas probed in absorp- 
tion. In this section, we also compare our results to modem theo- 
ries for star formation in atomic gas. Section 3 further extends our 
analysis beyond the inner few kpc where absorption arises, explor- 
ing the connection between DLAs and star-forming galaxies at any 
impact parameter. We also compare our findings with previous ob- 
servational studies. In Section 4, we use our findings to inform the 
discussion of what DLAs are, while a summary and conclusions 
follow in Section 5. 

Throughout this work, consistent with our previous analysis, 
we express distances in proper units and magnitudes in the AB 
system, adopting the following cosmological parameters: Ho = 
70.4 km Mpc“^, flm = 0.27 and Qa = 0.73 (Komatsu et 
al. 2011). As in previous papers of this series, we refer to the ab- 
sorption systems as “DLA gas” or simply DLAs, while we will 
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Table 1. Measured in-situ SFRs in different apertures centered at the quasar position for individual fields and for the composite images. 


Field 

^dla 

Dlum 

(Mpc) 

^igm 

(mag) 

/lin,(< 0.25”) 
(lO-^^erg s-lcm-^Hz-l) 

4>{< 0.25”) 
(Msyr-i) 

1.5”) 

(10“^^erg 

V>(< 1.5”) 
(M0yr“l) 

1:G1 

2.9181 

25152 

1.43 

- 

- 

<6.85 

< 1.05 

2:G2 

2.6878 

22760 

0.62 

- 

- 

<4.53 

<0.60 

3:G3 

2.3887 

19706 

0.18 

- 

- 

2.6± 1.0 

0.3± 0.1 

4:G4 

2.6878 

22760 

0.62 

- 

- 

<4.92 

<0.65 

5:G5 

2.5713 

21563 

0.20 

- 

- 

11. 7± 2.6 

1.4± 0.3 

6:G6 

3.7861 

34422 

1.80 

- 

- 

< 16.43 

< 3.85 

7:G7 

2.7544 

23449 

0.80 

- 

- 

< 4.70 

<0.65 

8:G9 

2.7584 

23490 

0.81 

- 

- 

< 3.49 

<0.49 

9:G10 

2.4592 

20420 

0.23 

- 

- 

< 1.84 

< 0.21 

10:G11* 

3.5297 

31646 

0.78 

- 

- 

7.2± 2.3 

1.5L0.5 

11:G12 

2.6606 

22480 

0.56 

- 

- 

< 3.93 

<0.51 

12:G13 

2.5978 

21835 

0.43 

- 

- 

< 3.29 

< 0.41 

13:H1 

2.7803 

23717 

0.23 

< 3.81 

< 0.54 

- 

- 

14:H2 

1.9127 

14990 

0.78 

<3.25 

< 0.24 

- 

- 

15:H3 

3.2530 

28685 

1.08 

< 2.68 

< 0.49 

- 

- 

16:H4 

2.7067 

22955 

0.20 

<3.76 

<0.51 

- 

- 

17:H5 

3.0010 

26021 

0.32 

<4.13 

< 0.66 

- 

- 

18:H6 

2.9586 

25576 

0.31 

<3.04 

< 0.48 

- 

- 

19:H7 

3.3069 

29259 

0.40 

< 1.59 

< 0.30 

- 

- 

20:H8 

2.6320 

22186 

0.31 

< 2.28 

<0.29 

- 

- 

21:H9 

1.8639 

14518 

0.58 

< 2.58 

< 0.18 

- 

- 

22:H10 

2.6826 

22707 

0.56 

<3.36 

< 0.44 

- 

- 

23:H11 

3.2332 

28474 

1.03 

< 1.88 

< 0.34 

- 

- 

24:H12 

3.5629 

32004 

0.88 

<2.99 

< 0.64 

- 

- 

25:H13 

2.6500 

22371 

0.39 

< 2.86 

< 0.37 

- 

- 

26:H14 

2.7333 

23230 

0.21 

<2.79 

<0.38 

- 

- 

27:H15 

2.9236 

25210 

0.30 

< 2.95 

< 0.45 

- 

- 

28:H16 

3.4035 

30291 

0.52 

< 1.81 

<0.36 

- 

- 

29:H17 

2.6289 

22154 

0.30 

< 2.04 

< 0.26 

- 

- 

30:H18 

2.6079 

21938 

0.25 

< 2.34 

< 0.30 

- 

- 

31:H19 

2.8721 

24672 

0.41 

<0.85 

<0.13 

- 

- 

32:H20 

2.6722 

22600 

1.02 

< 3.21 

< 0.42 

- 

- 

Stack 

- 

- 

- 

- 

< 0.090 

- 

< 0.270 


* For this sightline, we cannot exclude contamination from quasar light (see Section 2.1.1 for details). The columns of the table are (1) name of the quasar 
field; (2) redshift of the DLA; (3) luminosity distance to the DLA; (4) correction for intervening IGM absorption; (5,7) limit on the FUV flux measured at the 
DLA position, within the aperture specified in parenfhesis; (6,8) limit on the in situ SFR measured at the DLA position, within the aperture specified in 
parenfhesis. 


explicitly refer to the host galaxies as “DLA galaxy”. Furthermore, 
we will often refer to two sub-samples: the first one, the “HST sam- 
ple”, includes 20 DLAs which have been observed with WFC3 on 
board of FIST; the second one, the “ground-based sample”, includes 
12 DLAs which have been observed from the ground at Keck or at 
the Large Binocular Telescope. Details on the properties of these 
two sub-samples can be found in Fumagalli et al. (2014b). 


2 THE IN-SITU STAR FORMATION OF DLAs 

In this section, we take advantage of the unique possibility offered 
by the adopted double-DLA technique to explore, for the first time 
with direct measurements in the rest-frame UV, the star-formation 
activity associated with the absorbing DLA gas. We start by ana- 
lyzing emission properties in-situ, that is at the position of the ab- 
sorbing gas, inside individual DLAs (Section 2.1.1). We then gener- 
ate composite images to extend our study to fainter limits (Section 
2.1.2). Next, we connect these direct measurements to the absorp- 
tion properties, discussing in turn the implications for the star for- 


mation law (Section 2.2), the metal enrichment (Section 2.3), and 
the cooling rates (Section 2.4) of DLAs. 

Before proceeding, we emphasize that our observations probe 
FUV emission; thus, the measured SFRs represent only the fraction 
of the intrinsic star formation that originates in unobscured regions. 
In principle, we cannot rule out the existence of obscured star for- 
mation, but reddening measurements against quasars hosting DLAs 
imply the presence of only a modest amount of dust, with mean 
dust-to-gas ratios Ay /Nm ~ 2 — 4 x 10“^^ mag cm“^ (Vladilo 
et al. 2008; Khare et al. 2012). In the following, we will quote ob- 
served fluxes, non corrected for dust extinction. To help gauge the 
effects of dust extinction in our sample, we note that dust obscura- 
tion at the DLA position should lie in the range Av ~ 0.004 — 0.06 
mag in the DLA rest- frame. For an SMC-type extinction curve, this 
translates into a flux correction between 2—40% at 1200 A, depend- 
ing on the gas column density. This calculation reassures us that, at 
least to first order, the in-situ SFR measurements of this section are 
representative of the bulk of the ongoing star formation in DLAs. 
However, as we will discuss in more detail below, our measure- 
ments average emission from regions that are more extended than 
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the typical scales probed by absorption. Therefore, we cannot fully 
exclude the presence of clumps of obscured star formation which 
may lie in close proximity to the quasar sightline. 

2.1 The star formation rate in DLA gas 

2.1.1 SFRs in individual systems 

We start by integrating the flux within two apertures of diameters 
0 = 0.25” and 8 = 1.5” in each image. These apertures represent 
physical sizes of ~ 2 and ~ 12 kpc at the median redshift of the 
DLA sample^. These two apertures probe a population of DLAs 
that form stars in situ, either within compact or extended regions. 
We emphasize that in our analysis we consider the background 
quasar position as a reference, although the location of a sightline 
may not coincide with the geometric center of the physical struc- 
ture that gives rise to the DLA. The sky level is computed locally in 
annuli centered at the DLA positions, extending between 4” — 14” 
in radius. For the ground-based imaging, we compute SFRs only 
within the large aperture, which is more extended than the size of 
a point source. Conversely, we refrain from computing limits in 
the larger aperture for the HST sample, because, at the depth of the 
HST images, these limits are not informative compared to those ob- 
tained from the ground-based sample. These integrated fluxes, /int , 
are then compared to the flux upper limits, /ui, at 2cr confidence 
level (C.L.), which we compute by integrating the error models dis- 
cussed in Fumagalli et al. (2014b) inside the same apertures. The 
final flux limits are set to be fum = max(/int, /ui). Lacking spec- 
troscopic information, this definition allows us to express the most 
conservative upper limit on the in-situ SFRs of DLAs by treating 
positive detections as potential interlopers along the line of sight. 

As discussed in Fumagalli et al. (2014b), our sightlines are se- 
lected to ensure that the higher redshift LLSs (the natural blocking 
filters) fully block the quasar light at the wavelengths of our obser- 
vations. In fact, in most cases (20/32) the higher redshift absorber 
is a DLA or a super LLS, and in all cases it is a highly optically- 
thick (r > 2) LLS. For a typical quasar spectral energy distribution 
(SED) Fqso oc and given the redshift distributions of the in- 

tervening LLSs, the transmitted flux integrated in the filters used 
for this imaging survey is <?( 1% of the intrinsic flux for blocking 
filters with Nni = 10^® cm“^, and < 2% for a blocking filter with 
A^hi = 10^® cm“^. We therefore exclude leakage from the back- 
ground quasars as an important contaminant in our measurements. 
This will also become evident in the following discussion, given 
that only 3 DLAs show emission at the position of the absorber. Of 
these, two cases are clearly associated to a neighboring galaxy. 

The 10:G11 sightline deserves a particular note. From the 
gallery presented in Figure 8, one can see emission right at the 
quasar position. For this DLA, a strong LLS at ^us = 4.4671 acts 
as the blocking filter, but we can only constrain the column density 
to a lower limit of Nm > 10^^'® cm“^. Damping wings are not 
readily visible in the spectrum, but we cannot trivially set a robust 
upper limit, as this is a proximate system. For the observed quasar 
flux density Fqso ~ 4.7 x 10“^®erg at ~ 5000A, 

the upper limit of the quasar flux leakage within the imaging filter 
is Fqso < 1.7 X 10“®® erg As shown in Table 1, 

the flux detected against this quasar is below this limit, and there- 
fore we cannot exclude leakage as a source of contamination for 
this sightline. 

^ For the adopted cosmology, 1" corresponds to 8.1 kpc at z — 2.7. 


Throughout this work, fluxes (/i^) are converted into lumi- 
nosity (L^) and SFR (tp) with the calibration (cf. Fumagalli et al. 
2010b) 

Ip = 7.91 X 10“®®47rDium(^;)/i^Afcorr(^)(/igm(i;) (1) 

where Dium is the luminosity distance to the DLA, Kcon is the 
K-correction, and (pigm accounts for the fact that our observations 
probe flux blueward of 12 15 A in the DLA rest- frame and thus 
should be corrected for absorption in the IGM. We emphasize that 
the quoted SFRs bear significant uncertainty, especially given that 
we probe the FUV SED at < 1200A in the DLA rest-frame. At 
these wavelengths, UV fluxes are not an optimal tracer of recent 
star formation due to both a flattening of the SED around 1000- 
llOOA, and the onset of strong absorption lines (e.g. Leitherer et 
al. 2002). Nevertheless, to first approximation, a simple power law 
/x oc A“^ with /3 = 2 can be assumed to model the EUV con- 
tinuum between ~ 1000 — 1500A to within a factor of two. This 
error, albeit large, is comparable to the scatter in the EUV slopes 
from galaxy to galaxy and it is also comparable to the scatter in 
the model SEDs for different input metallicities and star formation 
histories. A power-law index 0 = 2 further allows us to adopt stan- 
dard calibrations at 1500A, and to trivially compute K-corrections. 

We estimate the IGM correction for each DLA by combin- 
ing the appropriate filter transmission curve with the mean IGM 
transmission, computed following standard procedures (e.g. Madau 
1995), but using the updated calculations by Inoue et al. (2014). 
We also account for the fact that the UV flux is fully absorbed at 
wavelengths blueward of the DLA Lyman limit. Additionally, ab- 
sorption from molecular gas and metals are considered negligible 
(Noterdaeme et al. 2008; Becker et al. 2013; Jorgenson et al. 2014). 
As discussed above, we do not correct for intrinsic dust absorption, 
which we argued is modest in DLAs. The final limits on the SERs 
are listed in Table 1 where formal 2cr upper limits are marked by 
the ’<’ sign to distinguish them from the few instances in which 
positive flux is detected. The listed errors account only for the un- 
certainty in the flux, without systematic errors on the SFR calibra- 
tion. 

The measured upper limits on the in-situ SFRs of individ- 
ual DLAs are shown in Figure 1, both for the HST sample (blue 
squares) and for the ground-based sample (red circles). Assuming 
the median upper limits as a typical reference value (dashed lines 
in the figure), we conclude from the HST imaging that DLAs do 
not form stars at rates >0.38 Mq yr“^ within compact regions of 
~ 2 kpc in size. Similarly, the ground-based sample indicates that 
typical DLAs do not form stars in extended regions of ~ 12 kpc 
in size with SFRs > 0.65 Mq yr“^. However, differently from the 
HST sample in which we do not have positive detections, for the 
ground-based sample we detect flux at more than 2 <t C.L. in 3/12 
apertures. As noted above, we cannot exclude that the flux detected 
against the DLA lOiGll is free from quasar contamination, thus 
conservatively we conclude that ^ 9% of the DLAs form stars with 
a rate above our sensitivity limits, and that virtually no DLAs form 
in-situ stars at these rates within compact clumps. Given the small 
sample size, these percentages bear significant uncertainty, but they 
offer the first quantitative description of the local star forming prop- 
erties of DLA gas in the high redshift Universe from FUV tracers. 

As for the positive detections, without spectroscopic redshifts, 
we do not know if these sources are genuine DLA galaxies. How- 
ever, for a given magnitude m! and projected distance r' from 
the DLA, we can compute the expected number N' of interlopers 
with m ^ m' and r ^ r' by using the observed galaxy num- 
ber counts in deep U and B band images. Using the published 
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Redshift 

Figure 1. Upper limits on the in-situ SFRs of DLAs in two apeitures of 
2 and 12 kpc (blue squares and red circles), computed respectively for the 
HST sample and for the ground-based sample. Limits are at 2fr for non- 
detections, and correspond to the measured flux for the three detections (see 
text for details). Filled symbols are for individual DLAs, while the dashed 
lines indicate the median upper limits of the two distributions. Upper limits 
measured in the composite images are shown instead with thick blue and 
red lines. 

values by Kashikawa et al. (2004) and Grazian et al. (2009), we 
find N' ~ 0.09, N' ~ 0.01, and N' ~ 0.0002 for 3:G3, 5:G5, 
and 10:G11 respectively. This translates into a low probability of 
finding random projected galaxies at the separations of 5:G5 and 
10:G11. Conversely, the probability that 3:G3 is an interloper is 
not completely negligible. Keeping the caveat of possible leakage 
in mind for 10:G11, it is plausible that some of these detections, 
and especially 5:G5, are physically associated to the DLAs (or 
lower redshift absorbers). Spectroscopic follow-up of these candi- 
dates is now necessary, although this represents a non-trivial task 
given their faint magnitudes (> 25.5). 

2.1.2 SFRs in composite images 

To probe even fainter SFRs, we generate two composite images by 
separately combining the 20 HST images and the 1 2 ground-based 
images. These stacks are generated after converting each image to 
physical units by applying the SFR calibration discussed in Section 
2.1.1. We also map angular separations into physical distances. In 
this way, we remove the redshift dependence from field to field, 
which is shown in Figure 1. To create the composite images, we 
then redistribute the SFRs measured in individual pixels onto a grid 
of pixel size 0.32 kpc and 1.1 kpc, for the HST and ground-based 
samples respectively. This pixel size, close to the native image pixel 
size at the median redshift 2 = 2.7, ensures a fairly uniform map- 
ping of the original images in the final stack, with ~ 19 and ~ 11 
independent pixels entering each pixel in the final composite for 
the HST and ground-based samples, respectively. Given that we are 
not interested in accurate image reconstruction, we adopt a simple 
shift-and-add algorithm that preserves the total flux. 

When producing stacks, a choice has to be made on what 
statistics should be used in combining images. Given the small sam- 
ple size, especially for the ground-based sample, the mean stack 
shows a significant numbers of positive detections across the field, 
due to the large number of galaxies which lie in proximity to the 
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Figure 2. Top. Spatial distribution of the median SFRs in a box of 40 kpc 
centered at the DLA position for the 20 quasar fields imaged with HST 
WFC3 UVIS. This map has been smoothed with a kernel of 3 pixels for 
visualization purposes. The circle represents the 2 kpc aperture used for 
the analysis. Bottom. The black solid line shows the histogram of the SFRs 
measured in individual pixels with size of 0.32 kpc. The SFRs that enter this 
histogram have been multiplied by a factor of 10 to improve the visibility of 
the full width. The red dotted lines show instead the distributions of SFRs 
measured within 10000 random apeitures, the sizes of which have been 
matched to the aperture shown in the top panel. 

quasars in these deep images. The corresponding flux distribution 
as a function of distance from the DLAs encodes interesting in- 
formation regarding the clustering of sources near the targeted ab- 
sorbers. We will consider the above in Section 3. In this section, we 
are interested in characterizing the typical in situ star forming prop- 
erties of the DLA gas. We hence choose to use the median statistic 
as a better estimator of the SFRs of typical DLAs, given that it is 
less sensitive to the bright pixels associated with the few positive 
detections within the regions of interest. During our analysis, we 
empirically reconstruct the noise properties of the two stacks rather 
than propagating the error map from individual images. The two 
median images are shown in Figure 2 and Figure 3. As in the pre- 
vious analysis, we compute the SFRs within two apertures of 2 and 
12 kpc diameters to constrain the in-situ star formation in compact 
and extended regions using the HST and ground-based imaging, 
respectively. We also measure SFRs within a larger aperture of 4 
kpc in the HST image to test for the presence of more extended 
emission in this sub-sample. 
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Figure 3. Same as Figure 2, but for the median of the 12 fields imaged from 
the ground. Here we show a region of 80 kpc centered at the DLA position. 
The image has been smoothed with a kernel of 2 pixel. Each pixel in the 
composite image has a size of 1.1 kpc. 


A proper assessment of the eiTor is critical to establish whether 
a low flux level is present in the stacked images. The pixel stan- 
dard deviation of the composite images, shown in the bottom panel 
of Figure 2 and Figure 3, is well-behaved, implying that our me- 
dian images are free from large-scale gradients. Furthermore, this 
pixel standard deviation, computed empirically across the image, 
is consistent with the typical pixel variance recovered from a boot- 
strap technique with 500 iterations. Therefore, we can use the pixel 
variance as our estimator of the noise in the composite image. To 
further account for correlated noise, we estimate the flux standard 
deviation for the science apertures by measuring the standard devi- 
ation of the recovered SFRs in 10000 apertures which we randomly 
position inside the composite images far from the quasar location. 

For the FIST sample, the distribution of SFRs in these 
randomly-distributed apertures is well approximated by a Gaus- 
sian with standard deviation 0.045 M© yr“^. Therefore, the SFR 
Ip = 0.011 Mq yr“^, measured in the 2 kpc aperture, is consistent 
with the range of SFRs detected in random apertures. We hence 
place a 2cr upper limit of t/> < O.OOOMq yr“^ on the median in 
situ SFR of the DLAs of the HST sample. For the ground-based 
sample, instead, the distribution of SFRs measured in random aper- 
tures is more skewed towards positive values. This is because, dif- 
ferently from the HST sample, both the smaller sample size and the 


higher density of detected sources bias the median of each pixel to- 
wards positive values (in a correlated fashion). Nevertheless, this 
distribution can still be approximated by a Gaussian with stan- 
dard deviation 0.14 Mq yr“^. It follows that the measured SFR 
tp = 0.13 Mq yr“^ within a 12 kpc aperture at the quasar position 
is statistically consistent with random noise, and we set a 2a up- 
per limit of Ip < 0.27 Mq yr“^ for the DLAs of the ground-based 
sample. These limits are summarized in Table 1 and Figure 1 . 

Finally, no flux is detected in the larger 4 kpc aperture in 
the HST composite, and we obtain a 2a upper limit of tp < 
0.23 Mq yr“^. Low surface brightness flux from even more ex- 
tended sources would be buried under the noise, which increases as 
a function of the number of pixels within the aperture. For compari- 
son, Reddy & Steidel (2009) report M* = — 20.97±0.14 at 1700A 
as the characteristic luminosity for star-forming galaxies between 
2.7 ^ z < 3.4, and thus our limits at ^ llOOA translate into 
SFRs which are typical for galaxies with ^ O.OIL, for the HST 
composite and ^ 0.03L* for the ground-based composite. 

Before proceeding further, we validate the adopted procedure 
by combining the HST images after artificially inserting at the DLA 
position in each image a source with size 1.6 kpc and SFR that is 
distributed as a Gaussian with median tp = 0.11 Mq yr“^. After 
processing these fake images through the same analysis pipeline, 
we recover an SFR t/> = 0.119 ± 0.045 Mq yr“^ in the composite 
image. Similarly, when we insert a source with size 5.5 kpc and 
SFR distributed around a median ip = 0.25 Mq yr“^ in each image 
of the ground-based sample, we recover a median ip = 0.38 ± 
0.14 Mq yr“^ from the stack. Thus, in both cases we recover the 
input SFRs within the errors. 

2.1.3 Comparison with previous work 

Both the analysis of individual systems and the study of the two 
composite images corroborate a picture in which DLAs do not form 
stars in situ with rates ip > 0.5 Mq yr“^. With the higher resolu- 
tion HST images, we also exclude the presence of compact star 
forming regions with SFRs as low as ip ~ 0.1 Mq yr“^ at the po- 
sition of the quasar. Keeping in mind that the conversion between 
the observed UV fluxes and the SFRs are subject to a substantial 
degree of uncertainty, we can put these limits in the context of pre- 
vious attempts to establish the SFRs of DLAs. 

By stacking spectra of DLAs with logA^ni ^ 20.62 and 
examining the residual flux in the Lyo trough, Rahmani et al. 
(2010) set a 2a upper limit on the Lya emission from DLA 
gas at < 2.0 x 10“^® erg cm“^, which corresponds to 
< 0.02L, at z ~ 2.8 for their assumed Lya luminosity func- 
tion, orip< 0.8 Mq yr“^ given their adopted conversion between 
Lya flux and SFR. From a re-analysis of the same data, Rauch & 
Haehnelt (2011) concluded instead that flux is positively detected 
to (5.35±1.97) X 10“^® erg cm“^ when restricting the analy- 
sis to positive velocities compared to the systemic DLA redshift, as 
expected for resonant lines. However, the detection of Lya in this 
composite spectrum is complicated by subtle systematic errors and 
it is only marginal. 

In fact, more recent work by Cai et al. (2014) sets an even 
more stringent limit on the Lya emission from a composite of 
~2,000 DLAs to < O.OIL* at z: ~ 2.6. Emission is detected when 
restricting to the highest column density DLAs with log A^ei > 
5 X 10^^ cm“^ (Noterdaeme et al. 2014), which are however a 
small subset of the general DLA population. Considering the many 
uncertainties at play in the two different measurements, we con- 
clude that both the analysis of Lya in composite spectra and our 
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Figure 4. Limits on the SFR surface densities Sgfr of DLAs from the HST 
sample (blue) and the ground-based sample (red). Limits for individual sys- 
tems are shown as squares and circles for the two samples, respectively. 
Limits derived from the composite images are shown as horizontal bars. 
Also shown with a grey shaded region is the range of SFR surface densi- 
ties detected in the outskirts of LBGs by Rafelski et al. (2011) and, with 
green dashed bands, the limiting surface densities probed by Wolfe & Chen 
(2006) for two different kernel sizes (0.25^^ for the upper value, and 1.0^^ 
for the lower value). The dotted black line shows where the extrapolation of 
the local SF law would lie at the lower column densities that are typical for 
DLAs. 


direct measurement of the FUV flux are in agreement and rule out 
the presence of appreciable in-situ star formation in DLA gas to 
comparable limits. 


2.2 The star formation law of DLAs 

Having established direct limits on the in-situ SFRs in a represen- 
tative sample of DLAs, we now turn to the relationship between 
the neutral gas observed in absorption and the rate with which 
new stars form, as probed in emission. Empirically, the link be- 
tween gas and star formation is described via a star formation (SF) 
law. One common parametrization of the SF law, known as the 
Kennicutt-Schmidt law (Schmidt 1959; Kennicutt 1998), is in the 
form Sstr = A'(Sgas/Sgas,o)^, where Egfr and Sgas are the SFR 
and gas surface densities, with constants K, Egas,o, and j3. While 
comparing Egn and Egas in our observations, we should bear in 
mind that pencil beam surveys intrinsically probe scales of few par- 
secs, while our SFR surface density extends to scales of few kpc. 
Therefore, this does not allow us to formally define in quantita- 
tive terms an SF law, which is known to be scale-dependent in the 
nearby Universe (e.g. Schruba et al. 2010). For this reason, and be- 
cause the measured upper limits do not allow us to constrain the 
SF law parameters in detail, here we only present a general discus- 
sion of the link between Egfr and Egas, also in comparison with 
previous studies. 

Figure 4 shows that, once recast as SFR surface densities, our 
limits become logEgfr < —0.90 Mq yr“^ kpc“^ in the compact 
regions probed by the HST imaging, and log Egfr < —2.25 Mq 
yr“^ kpc“^ in the more extended regions probed by the ground- 
based imaging. The composite images provide even tighter con- 
straints on the median SFR surface densities of the DLA popu- 
lation: logEsfr < —1.54 Mq yr“^ kpc“^ for the HST sample, 


and log Esfr < —2.62 Mq yr“^ kpc“^ for the ground-based sam- 
ple. Once compared to the Hi column densities, we see that DLAs 
could, in principle, consume the available neutral gas reservoir 
quite rapidly, despite their modest in-situ SFRs (cf. Prochaska & 
Wolfe 2009). 

For a range of DLA column densities that brackets our sam- 
ple (AThi = 10^°'^ — 10^^'^ cm“^), we can infer lower limits on 
the depletion time scales, which lie in the range > 0.06 — 0.56 
Gyr given the median SFR limits of the HST sample. Similarly, for 
the ground-based sample, the depletion time scales lie in the range 
> 0.66 — 6.66 Gyr. These values are simply a reflection of the 
fact that for typical DLAs, our limits lie for the most part above 
the extrapolation of the local SF law considering for example the 
original fit by Kennicutt (1998) with /3 = 1.40 ± 0.15, iLdisk = 
(2.5 ± 0.7) X lO"'' Mq yr“^ kpc“^, and Egas.o = 1 Mq pc“^. 
As shown by Figure 4, only the SFR limits for the highest column 
density DLAs in the ground-based sample are comparable to the 
extrapolation of the local SF law. Therefore, deeper limits, by 1 — 2 
orders of magnitude, would be required to probe the in-situ SFRs 
if the DLA gas had to form stars according to the extrapolation of 
the local SF law. This finding is in line with previous studies, as 
discussed in Section 2.2.1. 

Given the recent progress in understanding the phenomenol- 
ogy and the physics of the SF law, we should also consider from 
a theoretical point of view whether some star formation, if any, is 
to be expected in DLA gas with column densities Aini = — 

10^^ ® cm“^. Modem studies (e.g. Wong & Blitz 2002; Fumagalli 
& Gavazzi 2008; Bigiel et al. 2008; Fumagalli et al. 2010a; Bigiel et 
al. 2010; Gnedin & Kravtsov 2010; Ostriker et al. 2010; Feldmann 
et al. 2011; Krumholz et al. 2012; Glover & Clark 2012; Leroy 
et al. 2013) have shown that new stars preferentially form in re- 
gions of higher gas surface density and metallicity, where gas can 
cool efficiently. Conversely, star formation is reduced or completely 
suppressed in regions of low metallicity and at low column densi- 
ties. More quantitatively, both simple numerical models or detailed 
hydrodynamic calculations (e.g. Krumholz et al. 2009b; Gnedin & 
Kravtsov 2010) reveal the existence of a threshold in Egas.t below 
which Esfr sharply drops. At the low metallicity of DLA gas with 
Z ~ 0.1.^q, Egas,t ~ 10 Mq pc 

While the exact relationships between metallicity, molecular 
gas, and star formation are still being investigated (e.g. Glover & 
Clark 2012; Krumholz 2012), high in-situ SFRs in DLAs should 
not be common given that these systems lack a favorable envi- 
ronment for star formation (cf Krumholz et al. 2009a). Deeper 
imaging surveys that use the same strategy adopted here have the 
potential of directly constraining the physics of star formation in 
the lower column density and lower metallicity regions common to 
DLAs, offering a direct test of these models. This will be possible 
with future 30m telescopes, but it is also within the reach of current 
facilities, provided one carefully selects targets in narrow redshift 
ranges so as to maximize the FUV emission within the imaging 
filters. 


2.2.1 Comparison with previous work 

A few studies have already investigated the star-formation law in 
DLAs, using UV as a tracer for star formation (e.g. Wolfe & Chen 
2006; Rafelski et al. 201 1). However, because of the bright quasar 
glare, these analyses rely on indirect or statistical methods to link 
the Hi column density seen in absorption to the UV emission. The 
first of these studies is the work by Wolfe & Chen (2006), who 
searched for extended low-surface brightness galaxies with sizes 
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between 2 — 31 kpc in the Hubble Ultra Deep Field (HUDF). This 
search yielded for the most part non-detections, which were used 
by the authors to constrain the comoving SFR densities of DLAs 
at^ = 2.5 — 3.5. Comparing these limits to the SFR densities ex- 
pected if DLAs were to form stars according to an SF law, Wolfe & 
Chen (2006) concluded that star formation in DLAs proceeds with 
an efficiency^ of less than 10% compared to nearby galaxies. More 
recently, the study by Rafelski et al. (201 1) revisited the problem of 
the SF law in DLAs, under the assumption that these absorbers arise 
instead in the outskirts of compact LBGs. By stacking a sample of 
~ 50 LBGs at 2 ~ 3 in the HUDF, Rafelski et al. (2011) detected 
FUV emission to radii of ~ 6 — 8 kpc from the LBG centres. After 
connecting the observed SFR densities to the column density dis- 
tribution of DLAs, they concluded that the SFR efficiency of this 
gas is a factor of 10 — 50 lower than what one would predict from 
the local SF law. 

Figure 4 offers a comparison between our survey and these 
previous results, after rescaling all SFRs to the calibration assumed 
in this work. Our composite images probe SFR surface densities 
which are comparable to the ones probed by Wolfe & Chen (2006) 
for both compact and more extended regions. However, the orig- 
inal work of Wolfe & Chen (2006) does not provide direct limits 
on the SFR surface density of DLAs, but only limits on the SFR 
comoving density of the DLA population. Thus, our analysis com- 
plements and extends the results of their investigation, as it pro- 
vides direct observational evidence that in fact DLAs do not form 
stars in-situ within regions of appreciable UV surface brightness. If 
these limits are extrapolated to all DLAs up to a column density of 
AThi~ 10^^ cm“^, our observations naturally explain the lack (or 
the scarce number) of low-surface brightness sources in the HUDF, 
reinforcing the conclusions of the statistical study of Wolfe & Chen 
(2006). 

Compared instead to the analysis of Rafelski et al. (2011), lim- 
its from our composite HST image are not informative, given the 
superior depth of the HUDF and the larger sample size of LBGs 
used by these authors to create their stack. Our deep ground-based 
imaging, however, places more interesting constraints on the SFRs 
in an extended low-surface brightness component, similar to the 
one probed around LBGs by Rafelski et al. (2011). At first, one 
may conclude that the limit on the median SFR surface density 
which our study recovers is in tension with the surface brightness 
detected in the LBG composite by Rafelski et al. (2011). However, 
this is not the case for two reasons. 

First, Rafelski et al. (2011) produced a stack of compact 
LBGs, centering at the position of the peak surface brightness. 
Thus, for an underlying radial surface brightness profile, faint flux 
levels are coherently stacked when moving away from the LBG 
center. Conversely, in our stack, we coadd images centering at the 
DLA position, thus off-centre compared to the core of a putative 
nearby LBG. This has interesting consequences for the expected 
number of LBGs within a given projected distance from the quasar, 
and we will return to this point in the next section. For the pur- 
pose of this discussion, we note instead that, because DLAs can 
arise from many impact parameters, different surface brightness 
levels are not coadded coherently in our stack. This geometric ef- 
fect causes an intrinsic dilution of the signal, and thus our limits 


^ In agreement with previous DLA studies, in this section, we define the 
constant K in the star-formation law as efficiency, noting that modern in- 
vestigations on the SF law prefer the definition of depletion time for the 
inverse of this quantity. 
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Figure 5. Lower limits on the time that is needed for the m-sitii star forma- 
tion to enrich the DLA gas to the observed levels. As in previous figures, we 
show limits derived in 2 kpc apertures from HST imaging with blue squares, 
and limits derived in 12 kpc apertures from ground-based imaging with red 
circles. The green circle and the green square refer to the limits obtained us- 
ing the composite images in each of the two samples. The dark grey shaded 
region indicates the elapsed cosmic time between z = 10 and the DLA 
redshifts, which we assume as the maximum time available for the enrich- 
ment. This calculation relies on the assumption that the Ah i measured in 
absorption is representative of the Hi column density within the adopted 
apertures. 

are not necessarily at odds with the measurement of Rafelski et al. 

(2011). 

Secondly, and perhaps most importantly, according to the 
model proposed by Rafelski et al. (2011), the detected emission 
arises from column densities Nhi > 10^^ cm“^, which are not 
typical of the DLAs in our sample. For this reason, and given the 
above discussion about the SF law, the measurement of Rafelski et 
al. (201 1) applies to only a small fraction of the DLA population. 
Intriguingly, a similar trend has been seen in studies of Lya emis- 
sion in the DLA trough, where flux is detected in composite spectra 
of the highest and most rare column density DLAs (Noterdaeme et 
al. 2014), but not when considering the more typical lower col- 
umn density systems (Cai et al. 2014) included in our study. Our 
measurement and the previous discussion on the SF law thus high- 
lights how the results of Rafelski et al. (201 1) should not be simply 
extrapolated to the more general DLA population, outside of the 
column density interval that has been explored by our study. 

2.3 Implications for metal production 

By leveraging the new vantage point offered by our observations, 
we can compare directly the observed in-situ star formation of 
DLAs in emission to the observed metal content of DLAs in ab- 
sorption. During this part of the analysis, we restrict to the sub- 
sample of 26 DLAs for which we have high quality spectroscopy 
to compute the neutral gas phase metallicity, as discussed in Fu- 
magalli et al. (2014b). As emphasized for the previous analysis of 
the SF law, the absorption spectroscopy and the available imag- 
ing do not constrain the physical properties of DLAs on the same 
scales. For instance, the recent analysis by Kanekar et al. (2014) 
shows how the H I column densities measured on the smallest scales 
probed in absorption differ from the values measured when smooth- 
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ing on much larger spatial scales. Thus, the results presented in this 
section rely on the crude assumption that the measured column den- 
sities of hydrogen and metals are, to first order, representative of the 
typical DLA gas on scales of ~ 2 — 12 kpc. 

To link the limits on the SFRs and the measurements of metal- 
licity, we compare the mass in metals within DLAs (Mz.obs) to the 
upper limits on the mass in metals produced by in-situ star for- 
mation (Mz.sfr), given the constraints from our imaging observa- 
tions. With this comparison, we infer lower limits on the enrich- 
ment time Afz, a quantity which in turn offers a simple metric 
for whether the in-situ SFR alone is potentially sufficient to en- 
rich DLAs, or whether external contributions are needed. Specif- 
ically, we compute the mass in metals within DLAs as ALz.obs = 
lO^^/^^Z0(mpiVHi7rr^), where [X/H] is the metallicity measured 
in absorption and Zq = 0.0181 is the metallicity of the Sun (As- 
plund et al. 2009). In this equation, the term in parenthesis trans- 
lates the H I column density to the total hydrogen mass, assuming 
a cylinder of diameter r = 2 and r = 12 kpc for the HST and 
ground-based imaging, respectively. Given the negligible ioniza- 
tion of DLAs, the neutral gas phase is assumed to trace the bulk of 
the mass in the system. For this reason, our calculation is restricted 
to metals within DLAs, without considering metals in an ionized 
phase, which may nevertheless be physically associated to these 
absorbers (Fox et al. 2007; Lehner et al. 2014). 

Next, we compute the expected mass in metals related to 
star formation as Mz.sfr = yztpAtz, where yz is the metal 
yield, weighted according to a given initial mass function. Several 
choices for yz are available in the literature (see, e.g., the discus- 
sion in Peeples et al. 2014), with modem values ranging between 
~ 0.02 — 0.04. In the following, we assume yz = 1/42 (Madau 
et al. 1996), a common choice for such studies. Finally, comparing 
the inferred lower limits on Mz.sfr with the values of Mz.obs, we 
derive lower limits on the enrichment time, as shown in Figure 5. 

Flaving only upper limits on the SFRs, we cannot set strin- 
gent constraints on the history of the metal enrichment of DLAs. 
It is however interesting to note that the amount of metals locked 
in the DLAs is fairly modest once integrated over the 2 kpc area 
used for the analysis of the FIST sample. Thus, if the DLAs of our 
sample have SFRs just below the detection limit of the HST imag- 
ing (i.e. ~ 0.1 — 0.6 Mq yr“^), it would be straightforward to 
achieve the observed enrichment levels within ~ 0.05 — 10 Myr. 
For the ground-based sample, integrating over the larger (12 kpc) 
apertures, ~ 1 — 300 Myr would be needed to achieve the ob- 
served enrichment levels, if the SFRs were at the detection thresh- 
old (~ 0.2 - 2 Mq yr“^). 

Given that the available time for star formation to enrich DLA 
gas is ~ 2 Gyr if stars form as early as 2 : ~ 10, SFRs as low as 
Ip = 2 X 1O~^M0 yr“^, 500 times lower than the limits derived in 
our HST composite image, can still satisfy the observed metallicity 
in compact DLAs. A similar result holds for more extended DLAs, 
if star formation proceeds with a rate of ■0 = 9 x 1O~^M0 yr“^, 30 
times less than the limits of our ground-based composite image. We 
therefore conclude that the metals produced by low levels of in-situ 
star formation are in principle sufficient to aceount for the modest 
metal content of DLAs. However, if in fact no star formation oc- 
curs at the lower column densities common in DLAs, then external 
sources of metals are required to raise the mean DLA metallicity 
to ~ O.IZq, which is significantly above the metal content of the 
IGM (Schaye et al. 2003; Simcoe 2011). Again, deeper versions of 
this experiment have the potential of better constraining the mech- 
anisms through which the densest regions of the high-redshift Uni- 
verse are enriched. 


2.3.1 Comparison with previous work 

Previous work to constrain the enrichment histories of DLAs has 
compared two statistical quantities, i.e. the cosmic density of met- 
als produced by star forming galaxies (pz.sfr) and the amount of 
metals that reside in DLA gas at any given time (pz,dia). The con- 
sensus is that DLAs at z ~ 2 — 3 contain only a small fraction 
of the metals produced in star-forming galaxies (e.g. Pettini et al. 
1999; Bouche et al. 2007; Rafelski et al. 2014), with the latest es- 
timates suggesting pz.dia ~ O.Olpz.sfr at z ~ 2.3 (Rafelski et 
al. 2014). Albeit with substantial uncertainties, most notably not 
accounting for a potential contribution from an enriched and dust- 
obscured DLA population, it clearly appears that DLAs are not en- 
riched by all the metals produced in galaxies. A similar “missing 
metal problem” was found by Wolfe et al. (2003b), who compared 
pz.dia with a new estimate of pz.afr in DLAs (and not LBGs) de- 
rived from the cooling rates inferred from [Cll*] absorption lines 
(see also Section 2.4). 

However, the comparison between pz.dia and pz.sfr from in- 
direct constraints on the in-situ SFRs in DLAs has also yielded 
contrasting results, depending on the different assumptions on the 
nature of DLAs which were made by different authors. In particu- 
lar, combining their inferred star formation efficiency for extended 
low-surface brightness DLAs with the results of numerical simu- 
lations, Wolfe & Chen (2006) suggested that metals may be un- 
derproduced in DLAs compared to the observed values. When ac- 
counting instead for the yields associated with the extended UV 
emission detected in the outskirts of LBGs, Rafelski et al. (2011) 
achieved instead a closure of the metal budget in DLAs. 

The analysis presented in Figure 5 is consistent with the sta- 
tistical arguments at the origin of the missing metal problem. Given 
the lack of appreciable local star formation, as pointed out by our 
observations, one can understand why DLAs contain fewer met- 
als than those produced in LBGs, which form stars at much higher 
rates. Furthermore, as shown in Figure 5, DLAs can be enriched to 
the observed values even by low levels of star formation, either in- 
situ or, if tp 10“^ — 10“'^ M 0 yr“^, by modest episodes of star 
formation occumng in the immediate surroundings (e.g. Rafelski 
et al. 2011). Conversely, at the limits of our ohservations or those 
by Wolfe & Chen (2006), it seems unlikely that we are facing a 
problem of underproduction of metals in DLA gas. 


2.4 Implications for cooling rates 

The [Cll*]A1335.7 absorption line arises from the ^Ps /2 level in 
the ground state of ionized carbon (C*"), and thus provides a tool to 
measure the column density of C*" ions in the j = 3/2 state. This 
is of much interest as these ions give rise to the [CII] 158pm tran- 
sition (^Pa /2 — P 1 / 2 ), which is believed to be the major coolant 
of the neutral ISM (e.g. Wolfire et al. 1995). The strength of the 
[Cll*]A1335.7 absorption can hence be used to measure the cool- 
ing rate (ic) of DLA gas (Pottasch et al. 1979; Wolfe et al. 2003a), 
and, under the assumption of thermal balance, to infer the corre- 
sponding heating rate. Leveraging this idea, Wolfe et al. (2003a,b) 
have developed a formalism to quantify the intensity of the radi- 
ation field that is responsible for the photoelectric heating of the 
DLA gas, and hence for the corresponding SFR surface density. 

Applying this formalism to a sample of 38 DLAs for which 
Ssfr could be inferred, Wolfe et al. (2008) noted a bimodal dis- 
tribution in the inferred star formation surface densities, which 
reflects the bimodality in the DLA cooling rates. Measurements 
of £c show two different populations, one of “high-cool” sys- 
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Figure 6. A comparison between the SFR surface densities infeixed by 
Wolfe et al. (2008) from [C II*] absorption (histograms) and the limits in- 
ferred from our composite images (vertical lines). From left to right, the 
grey solid, green dashed, and brown dot-dashed histograms refer respec- 
tively to low cool DLAs, lower limits for the high cool DLAs, and values 
corrected for the filling factor of the star forming regions for the high cool 
DLAs. The median of the observed sample (low cool and limit on high cool) 
is also shown with a black triple dotted line. The Wolfe et al. (2008) values 
have been corrected by a factor 1.58 to account for the different SFR cali- 
bration used in our study. The vertical lines are for the SFR surface density 
limits measured in the 2 kpc aperture in the HST composite image (blue 
line to the right) and in the 12 kpc aperture in the ground-based composite 
image (red line to the left). 


terns and one of “low-cool” systems, separated at around Ic ~ 
erg H“^(see also Neeleman et al. 2013). While a low 
cooling rate can be balanced by the heating rate from the meta- 
galactic UV background alone, Wolfe et al. (2008) concluded in 
their study that an additional heating input is generally required 
by their model to describe these observations. These authors there- 
fore computed the SFR surface densities needed to satisfy the mea- 
sured cooling rates, under the assumptions detailed in Wolfe et al. 
(2003a,b). These distributions are shown in Figure 6 for both the 
low-cool (solid gray histogram) and the high-cool (green dashed 
histogram) DLAs. 

The ability to measure the UV surface brightness in the same 
DLAs with detected C It* absorption would offer a direct test of this 
model, in turn providing interesting constraints on the cooling and 
heating mechanisms at work in DLA gas. Unfortunately, our survey 
is designed to target 2 ~ 2.7 DLAs in 2 > 3.5 quasars, for which 
the C II* absorption is redshifted inside the Lya forest, hampering 
a reliable determination of the cooling rates in individual systems. 
Flowever, we can proceed with a statistical approach, by comparing 
in Figure 6 the limits on the in-situ SFR surface densities measured 
in the composite FIST and ground-based images, with the surface 
densities inferred by Wolfe et al. (2008) from the observed C ll* 
cooling rates. 

Considering the low-cool DLAs at first, we see that even the 
stringent limits derived from the composite images do not provide 
interesting constraints. A factor of 10 deeper survey is needed to 
probe the nature of this sub-population. Considering the high-cool 
DLAs instead, the null detection in our ground-based composite 
image appears at odds with the San implied by the high-cool DLA 
populations for extended galaxies. As shown in Figure 6, the me- 


dian Esfr in the ground-based sample is about 0.5 dex lower than 
the median of the Sgtr inferred from the cooling rates when includ- 
ing both low-cool and high-cool, under the assumption that our 
sample is not biased towards one of the two populations (but see 
below). In their work, Wolfe et al. (2008) reached a similar conclu- 
sion on how a population of star forming DLAs which satisfy the 
Ssfr required by the high-cool population would violate the SFR 
surface density limits imposed by the Wolfe & Chen (2006) anal- 
ysis. This is not surprisingly, given the previous discussion on the 
limits shown in Figure 4. 

As noted by Wolfe et al. (2008), this tension can however be 
alleviated if compact LBGs in proximity to, but not fully embedded 
with, the DLAs act as heating sources. Because of simple geomet- 
rical arguments related to the r~^ dependence of the flux, how- 
ever, the inferred SFR surface densities for the high-cool popula- 
tion should be corrected upward by the relative tilling factor of star 
forming regions and DLA gas. Once this correction is applied, the 
disti'ibution of Esfr for the high-cool DLAs shifts to higher values, 
as shown in Figure 6 with a brown dot-dashed histogram. Flow- 
ever, our observations constrain this possibility: if the heating aiises 
from compact LBGs embedded in DLA gas, these sources must be 
located at distances >1 — 6 kpc (i.e. the radii of the adopted aper- 
tures throughout our analysis) from the DLA position, as indicated 
by the lack of emission at Esfr O.IMq yr“^ kpc“^ in the FIST 
and ground-based samples (Figure 4). It is however possible that 
compact LBGs with high Esfr lie at lai'ger impact parameters, a hy- 
pothesis we will explieitly test for in Section 3.4 when studying the 
distribution of star-forming galaxies in proximity to DLAs. 

Before proceeding, we caution the reader that the basis for 
the comparison presented in this section is that, if the bimodal dis- 
tribution of the cooling rates and the inferred Esfr are typical of 
DLAs, than comparable SFRs should be detected in our sample 
which is representative of the DLA population. Flowever, it should 
be noted that the cooling rates of DLAs correlate with metallic- 
ity, with the majority of high-cool systems having [X/H] > —1.2 
(Wolfe et al. 2008). The distribution of metallicity of DLAs in our 
study peaks instead around [X/H] ~ —1.5 (see figure 14 in Fuma- 
galli et al. 2014b). Furthermore, throughout this analysis, we are 
assuming an equal split between the low-cool and high-cool pop- 
ulation, following Wolfe et al. (2008). However, the estimate of 
Wolfe et al. (2008) is based on a small sample of DLAs and may 
not reflect the actual distribution in the general DLA population. 
Our non-detection could therefore be simply explained by a larger 
percentage of low-cool DLAs, for which we do not expect UV flux 
in excess to our sensitivity levels. For these reasons, and also be- 
cause of the small number of systems in our sample, any tension 
between our observations and the inferred Esfr from [C ll*] absorp- 
tion should be regarded only as tentative at this point. Future work 
in larger samples, and in particular a direct comparison of the ob- 
served versus inferred SFR surface densities in individual systems, 
is needed to confirm this appai'ent tension. 


3 THE CONNECTION BETWEEN DLAs AND GALAXIES 

In Section 2, our discussion focused entirely on the emission prop- 
erties of DLA gas, directly at the position where spectroscopy in 
absorption reveals the presence of neutral hydrogen with column 
densities above Nm ^ cm“^. In this section, we expand 

our study beyond the properties of the in-situ SFRs of DLA gas, by 
considering associations between DLAs and star-forming galaxies 
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at all impact parameters 6dia, which we define here as the projected 
quasar-galaxy separation. 

Given that our images probe only the rest-frame FUV emis- 
sion, our analysis will only consider galaxies with appreciable un- 
obscured SFRs, typically ip > 1 — 2 Mq yr“^. Following early 
studies that employed color selection techniques to identify a pop- 
ulation of star-forming galaxies at z > 2.5 (e.g. Steidel et al. 1996; 
Madau et al. 1996; Giavalisco 2002), the term LEG has become 
a synonym for star-forming galaxy in the literature. Throughout 
this work, we follow the same convention. We clarify however that 
the term LEG here refers to any star-forming galaxy, regardless 
of its actual SFR. Traditionally, LEG samples included primarily 
UV-bright galaxies with SFRs of tens to hundreds of Mq yr“^, 
but our study also probes galaxies with lower SFRs, by a factor of 
10 to 100. In this discussion we will make a distinction between 
“bright” LEGs and “faint” LEGs, or dwarfs. With the expression 
bright LEGs, we will refer to a population with properties similar 
to the original samples of spectroscopically-confirmed LEGs, with 
SFRs t/) ~ 20 — 50 Mq yr”*^ and halo masses ~ 10^^ Mq. Con- 
versely, with the expression faint LEGs, we will refer instead to a 
population of star-forming galaxies with SFRs ip < 10 Mq yr“^, 
down to the sensitivity limits of our survey. For this second popu- 
lation, we will also distinguish between faint LEGs which are iso- 
lated, i.e. central to their dark matter halos, from faint LEGs which 
are satellites of brighter LEGs. The reason for this distinction will 
become apparent in the following discussion. 

Finally, we note that in the remainder of this paper, we will be 
naturally biased against a population of star-forming galaxies that 
is significantly dust obscured (see the discussion in Section 4.2). 
In Section 2, we argued that the low dust extinction seen along the 
quasar line of sight rules out severe dust obscuration as the cause 
of the observed lack of appreciable in-situ star formation. However, 
when considering galaxies at all impact parameters, we do not have 
any prior knowledge of the dust properties. Thus, the conclusions 
drawn from our following analysis apply to a galaxy population 
with unobscured SFRs comparable to our detection limits. Galaxies 
below this detection limit could in principle be still forming stars 
to appreciable intrinsic rates. Unfortunately, with our current data, 
we cannot assess the importance of dust obscuration in searches of 
DLA galaxies, which remains an open question that future studies 
at infrared and millimetre wavelengths should address. 

3.1 Limits on the impact parameters and SFRs 

Without a deep and dense redshift survey of the galaxies in the tar- 
geted quasar fields, we cannot establish unique DLA-galaxy asso- 
ciations for all the absorbers. However, we continue in the spirit of 
the previous analysis by studying the distribution of SFRs and im- 
pact parameters of the identified galaxies to derive statistical con- 
straints on the properties of the DLA galaxy population. We start 
our analysis by focusing on the galaxies in proximity to the quasar 
sightline (6 < 10 — 30 kpc). From a physical point of view, asso- 
ciations between DLAs and galaxies at these separations are use- 
ful to constrain models in which DLAs arise from a gaseous disk 
or a more amorphous structure in the surroundings of an embed- 
ded star-forming disk (e.g. Fynbo et al. 2008; Rafelski et al. 2011; 
Danovich et al. 2014). More practically, this is the interval of im- 
pact parameters for which our observations become the most con- 
straining, given the low number of random interlopers predicted 
within the small projected area subtended by impact parameters of 
6 < 30 kpc. 

To derive limits on the SFRs and impact parameters of DLA 


galaxies, we produce for each quasar field k a ranked list of im- 
pact parameters with associated SFRs {fedia.i, ipi)k, assuming that 
each galaxy i is at the DLA redshift. For this analysis, we limit our 
search to bdia = 200 kpc, which coiresponds to ~ 2 times the 
virial radius of a massive LEG at ^ ~ 2.7 (e.g. Fumagalli et al. 
2014a). This distance is considered sufficient to encompass both 
DLAs that arise from the immediate surroundings of galaxies and 
systems that originate from gas structures (including but not limited 
to tidal debris or accreting filaments) associated with the halo of a 
galaxy. In this ranked list, the impact parameters can be generally 
considered as lower limits of the true impact parameters because, 
if the i— th galaxy is not the DLA host, the i + 1 galaxy will be- 
come the next candidate with 6dia,i-n > 6dia,i- However, if the 
i— th galaxy is in fact the true DLA galaxy, or the DLA galaxy is 
fainter than our detection limit, then 6dia,j for j > i (or even for 
each j) will be distributed as dictated by random galaxy counts. It 
follows that, by comparing the observed {6dia,i, distribution 
in the quasar surroundings with the distribution of random galaxies 
{6ran,m, we Can identify which galaxies have impact param- 

eters and SFRs which are unusual given random galaxy counts, thus 
constraining the impact parameters and SFRs of the DLA galaxy 
population. 

We present results from this comparison in Figure 7, where 
we show the distribution of impact parameters and SFRs for all the 
galaxies detected as a function of the distance from the 32 DLAs 
included in this study. In each quasar field, we also select three 
apertures of 200 kpc in radius, which we place at random positions 
far from the quasar locations. All the galaxies detected in these ran- 
dom apertures are used to characterize {bran.m, 'tpra}k, that is the 
SFR and impact parameter distributions of a “background” galaxy 
population given random counts of galaxies that are not physically 
associated with the targeted DLAs. These distributions, which are 
shown in the right panels of Figure 7, are subject to the same se- 
lection effects as the galaxy population near the DLAs (see a dis- 
cussion in Fumagalli et al. 2014b). Thus, by empirically measuring 
the background in proximity to the regions of interest with the same 
selection function, we automatically account for the field-by-field 
variation in the sensitivity limits of our survey, the differences aris- 
ing from the fact that we have adopted multiple filters, and potential 
differences in the large-scale distribution of galaxies in the targeted 
fields. 

The main result of this analysis is already apparent from the 
histograms in Figure 7. The distributions of SFRs and impact pa- 
rameters for the galaxies close to the quasar sightline match the cor- 
responding distribution for random galaxies far from the quasars, 
without an obvious excess that can be attributed to the DLA hosts. 
The only noticeable difference is at 6dia < 10 kpc, where few 
galaxies are detected. The same excess is visible also as a func- 
tion of the SFR, as shown in the left panel of Figure 7. Even with 
these few detections, consistent with our previous discussion on 
the in-situ SFR, the region feaia < 10 kpc is almost completely 
empty, with only 3/32 detected galaxies (2/32 if we interpret the 
detection in the 10:G1 1 field as due to quasar leakage). From this 
analysis, we draw an important conclusion: not only do DLAs not 
form appreciable stars in situ, but also that < 7% of the DLAs are 
associated with galaxies with ip ^ 2 Mq yr“^ and 6dia < 10 kpc. 
In turn, this implies that the majority of DLAs arise from galaxies 
with Ip < 2 Mq yr”’^ and/or at 6dia > 10 kpc. 

Figure 7 also offers a convenient empirical way to summa- 
rize the effective completeness of our survey. One can see that 
the distribution of impact parameters reaches its maximum density 
around ~ 2 Mq yr“^, while it is more sparse at both larger and 
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Figure 7. Left: Distribution of impact parameters as a function of SFRs for the galaxies near DLAs. The grey crosses are for galaxies detected near the DLAs, 
and values are computed at the DLA redshift. Errors on SFRs are not shown to improve visibility, but they range between 5 — 20% (statistical). The red solid 
curve represents the median minimum impact parameters computed in 10000 trials by drawing from a “background” galaxy population, i.e. random galaxies 
detected far from the quasar position. The corresponding percentiles of these minimum impact parameters are shown with dotted and dash-dotted curves, 
respectively for the 10 — 90% and 25 — 75% of the distribution. The black curve shows a numerical approximation for the 25% of the distribution, which 
holds up to the completeness limit indicated by a vertical dashed line. Outlying galaxies (i.e. those lying below the expected minimum impact parameters for 
a random galaxy population) are labeled with the name of their quasar. Right: Normalized histograms of the SFRs (top) and impact parameters (bottom) for 
the galaxy population around the DLAs (black dashed histograms) and for the background galaxy population. 


smaller SFRs. At larger SFRs, such a decrease is expected accord- 
ing to the general form of the luminosity function, for which the 
number of sources decreases as the apparent magnitude decreases 
(hence the luminosity increases). The decrease of density for in- 
creasing apparent magnitudes is instead a result of incompleteness 
in our flux-limited survey. In the following analysis, we will con- 
sider ip = 2 M 0 yr“^ as our effective completeness limit, keeping 
in mind that the depth of individual images varies from field to 
field. We note that the discrepancy between the completeness limit 
and the more sensitive limits on the in situ SFR arises from the dif- 
ferent nature of the two measurements. While the in situ measure- 
ments are only limited by pixel variance, the completeness limit 
also depends on the procedures used to identify sources in the im- 
ages. We refer the reader to Fumagalli et al. (2014b) for a detailed 
description of these effects. 

Beyond 6 dia = 10 kpc, one can instead see an increasingly 
high number of detections, which, as suggested by the distribu- 
tions of SFRs and impact parameters, are for the most part unre- 
lated to the DLAs. But are we seeing any excess compared to a 
pure background distribution, which we can associate with the in- 
tervening DLAs? To answer this question more quantitatively, we 
derive an estimate for the locus of the minimum impact parameter 
as a function of SFR, given a random galaxy population and a ref- 
erence point in the sky. In other words, we seek for an expression 
6 min = hmin(t/’) which quantifies the typical distance to the closest 
galaxy in the sky, given an SFR and a preferred view point as set 
by a quasar. We emphasize that in this calculation ip is a surrogate 
for the apparent magnitude, and it should not be inteipreted as an 
actual SFR measurement. 


To derive &min, we simulate the locus of impact parame- 
ters as a function of SFRs using the distributions of background 
galaxies with 10000 mock samples. Specifically, for each trial, 
we draw iVran pairs {bra.n,b,ipb}k from the distributions of im- 
pact parameters and SFRs of the background galaxy popula- 
tion, with Aran equal to the number of sources detected within 
200 kpc from the 32 DLAs included in this study. We then 
define 6 min(V') = miii(^'ran,b(t/’)) ill bins of SFR. By re- 
peating this procedure 10000 times, we recover a distribution 
{&min}t, the median of which is shown in Figure 7, together 
with the 0.1,0.25,0.75,0.9 percentiles. To first approximation, 
^'min, 25 %/kpc = {ip /M q yr“^)®'® -|- 6 . This parametrization 
should be generalized to other samples with caution, being related 
to the extremes of a distribution. 

With an estimate for the typical bmin in a random galaxy sam- 
ple, we have a metric to conclude that there is a statistically sig- 
nificant excess of galaxies at small impact parameters, with 4/32 
detections below &min,io%> and 6/32 detections below bmin,26%- 
Images of these 6 quasar fields are shown in Figure 8 . Thus, from 
this analysis, we can conservatively conclude that ^ 13% of the 
DLAs are associated with galaxies with ip ^ 2 M© yr“^ and im- 
pact parameters fedia ^ bmin, 25 %- In Figure 7, we also label the 
location of the galaxy in the field 3:G3 that gives rise to the pos- 
itive detection within the 12 kpc aperture centered on the DLA in 
the previous section. This galaxy lies close to the typical bmin for 
random galaxies, in line with our previous estimate of a ~ 10 % 
probability to detect a random interloper with comparable magni- 
tude in the search area defined by the projected impact parameter. 
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Figure 8. Images in the bluest available filters of the six quasar fields with candidate DLA galaxies at close impact parameters. Each panel shows a 30" X 30" 
region centered at the DLA position, which is marked by a red circle. The detected sources are marked with white ellipses, which represent the Ki'on apertures 
used for photometry. 


The association of the detected emission for this DLA is therefore 
tentative. 

Nevertheless, we believe that we have detected a handful of 
sources that are likely associated to the intervening DLAs along the 
targeted sightlines. We emphasize that, although the target fields 
were selected due to the presence of a DLA and we are hence 
“biased” towards the association of any detected galaxy with the 
DLA, we cannot formally rule out associations with lower-redshift 
absorption line systems (e.g. Mgll or CIV absorbers) which are 
present in these spectra. This is especially true given the statistical 
nature of this argument and the small number of positive detec- 
tions. Thus, the 6/32 identifications should be regarded as quite 
conservative upper limits. 

Finally, we note that the host galaxies of the higher redshift 
LLSs and DLAs which act as blocking filters should not be de- 
tectable below their Lyman limit, unless they are characterized by 
an unusually high escape fraction of ionizing photons. Thus, the 
presence of higher redshift blocking filters does not affect our es- 
timates. In passing, we point out that the search for dropout galax- 
ies in proximity to these quasars will offer interesting candidates 
for the association with the higher redshift DLAs and LLSs in 
these sightlines. At the same time, the search for DLA hosts at 
&dia > 10 kpc through a statistical analysis like the one here 
presented does not have to be restricted to sightlines with double 
DLAs, and it can be easily performed in larger samples. We defer 
these searches to future work. 
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Figure 9. The projected DLA-galaxy correlation function as measured 
around the 32 DLAs included in this study within bins of 40 kpc {except 
the innermost bin of 20 kpc). The measurement obtained with real data is 
shown with black squares, while blue circles and red stars show the pro- 
jected correlation function after injecting 16 and 32 artificial DLA hosts 
with random impact parameters in the range — 20 — 100 kpc. For 
these two cases, we show only the two relevant bins, as the other values re- 
main unchanged. The central values have also been offset by 5 and 10 kpc 
for visualization purposes. 


3.2 The galaxy-absorber correlation function 

In section 3.1, we focused our attention on the closest impact pa- 
rameters, in the range baia ^ 30 kpc. We now extend our analysis 
to larger impact parameters to investigate whether there is an ex- 
cess of galaxies near DLAs compared to a background galaxy pop- 
ulation. To address this question, we measure the projected DLA- 


galaxy correlation function, by comparing the observed number of 
galaxies at a given separation from the targeted DLAs to the ex- 
pected number from a random galaxy population (see, e.g., Cooke 
et al. 2006). Again, for this calculation, we consider a circular aper- 
ture of 200 kpc around each DLA. 
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Table 2. Properties of intervening DLAs with confirmed galaxy associations. 


Name 

-2qso 

^dla 

ba.s 

n 

bp 

(kpc) 

log Ah I 
(cm-=) 

SFR 
(Mq yr“ 

-1) 

[XAI] 

Ref. 

QO 139-0824= 

3.017 

2.677 

1.60 ±0.05 

13.0 ±0.4 

20.70 ±0.15 

- 

- 

-1.15 ±0.15 

[16] 

Q0338-0005 

3.068 

2.230 

0.49 ±0.12 

4.1 ± 1.0 

21.05 ±0.05 

- 

- 

-1.25 ±0.10 

[16] 

PKS0458-02 

2.286 

2.039 

0.31 ± 0.04 

2.6 ±0.3 

21.77 ±0.07 

>1.5 

Lyo: 

-1.19 ±0.10 

[3,6,14] 

PKS0528-250 

2.797 

2.811 

1.14 ±0.05 

9.1 ±0.4 

21.27 ±0.08 

4.2 

Lyo: 

-0.75 ±0.10 

[4,5,9,16] 

Q0918+16361 

3.07 

2.583 

1.98 ±0.01 

16.2 ±0.1 

20.96 ±0.05 

8±3 

Ha 

-0.12 ±0.05 

[11,12] 

Q0953+47= 

4.457 

3.404 

0.34 ±0.10 

2.6 ±0.8 

21.15 ±0.15 

- 

- 

-1.80 ±0.30 

[16] 

J1 135-0010 

2.89 

2.207 

0.10 ±0.01 

0.8 ±0.1 

22.10 ±0.05 

25 ±6 

Ha 

-1.10 ±0.08 

[13] 

Q2206-1958a 

2.559 

1.920 

0.99 ±0.05 

8.5 ±0.4 

20.65 ±0.07 

5.7 

UV 

-0.54 ±0.05 

[7,9] 

Q2206-1958b 

2.559 

1.920 

1.2 ±0.1 

10.3 ±0.9 

20.65 ± 0.07 

4.2 

UV 

-0.54 ±0.05 

[2,9] 

Q2222-0946 

2.926 

2.354 

0.8 ±0.1 

6.7 ±0.8 

20.65 ± 0.05 

9.5 ± 1 

Ha 

-0.46 ±0.07 

[8,10,14,15] 

HE2243-60 

3.01 

2.328 

2.80 ±0.20 

23.4 ± 1.7 

20.62 ±0.05 

18 ±2 

Ha 

-0.72 ±0.05 

[1,16,17] 


References: [1] Bouche et al. (2013); [2] Weatherley et al. (2005); [3] Mpller et al. (2004); [4] Mpller et al. (1993); [5] Moller & Wairen (1998); [6] Wolfe et 
al. (2005); [7] Mplleretal. (2002); [8] Peroux et al. (2012); [9] Christensen et al. (2014); [10] Fynbo et al. (2010); [11] Fynbo et al. (2011); [12] Fynbo et al. 
(2013); [13] Noterdaeme et al. (2012b); [14] Krogager et al. (2013); [15] Jorgenson & Wolfe (2014); [16] Krogager et al. (2012); [17] Bouche et al. (2012). 

Notes: ^In the same sightline, Fynbo et al. (2013) reported the detection of [OIII] emission line associated to a second DLA at 2 — 2.412. ^This DLA has 
been reported by Krogager et al. (2012), citing work in preparation by other authors. Errors on the impact paramaters are from Krogager et al. (2012). We do 
not include in this list the proximite DLA Q0151+048A (Fynbo et al. 1999) and the sub-DLA 2233.9+1318 from Djorgovski et al. (1996). 


Specifically, the number of random galaxies inside an annulus 
A6 at a given projected separation b from a DLA is 

Afran — 27r6A6/tran , (2) 

where /Tran = (2.566±0.046) x 10“'^ kpc~^ is the surface number 
density of galaxies at the depth of our imaging survey, which we 
compute using the previously defined apertures at large distances 
from the quasar positions. By comparing A^ran to the observed 
number of galaxies Nobs in a given annulus around the DLAs, we 
obtain a binned estimator of the projected DLA-galaxy correlation 
function, Wdg(b), defined by 

Wdg{b, Ab) = Nobs{b, Ab)/Nr,n{b, Ab)-1. (3) 

Errors are simply set according to Poisson statistics in the observed 
number counts. A similar analysis, but employing mock galaxy cat- 
alogues constructed in apertures of 200 kpc and the DD/RR estima- 
tor, yields the same result. Given our qualitative discussion, we do 
not resort to the more advanced techniques for measuring correla- 
tion functions and estimating errors, which exist in the literature 
(e.g. Landy & Szalay 1993). 

The measurements of Wdg{b) are shown in Figure 9. In agree- 
ment with the nearest neighbor statistics presented in Section 3.1, 
the projected DLA-galaxy correlation function exhibits an excess 
for the innermost bin, in the interval 0 — 20 kpc. This excess is 
however only marginal, especially since our simple errors prob- 
ably underestimate the true variance. Furthermore, this excess is 
severely dependent on the choice of the bin size, as is typical for 
correlation functions (e.g. Croft et al. 1997). Increasing the inner- 
most bin size to 30 kpc would cause the signal to disappear due 
to the rapid increase in the number of background sources, as one 
could also guess from Figure 7. 

Conversely, for bjia > 20 kpc, we do not see any excess of 
galaxies in proximity to the DLAs compared to background values. 
This is expected, since we are computing the projected correlation 
function without spectroscopic information, thus integrating over a 
huge volume. We can gauge the sensitivity of our measurement to 
DLA host galaxies by recomputing the correlation functions after 
injecting 32 or 16 fake DLA galaxies in the observed catalogues 
with impact parameters haia = 20 — 100 kpc. In this simple exper- 


iment, we assume that all galaxies in the observed catalogues are 
random interlopers, and that the fake sources are the only galaxies 
associated to the DLAs. Results from this calculation are shown 
in Figure 9. As expected, the presence of tens of DLA galaxies at 
these impact parameters boosts the clustering signal, although the 
detection remains marginal. A comparison between Wdg{b) with 
mock and real data suggests that a large fraction (S> 50%) of DLA 
galaxies with SFR 'tp ^ 2M0 yr“^ and haia = 20 — 100 kpc are 
not present in our data, thus implying that not all DLAs originate 
nearby to highly star-forming galaxies. 


3.3 Comparisons with other DLA searches 

The search for DLA host galaxies has been a long-time effort 
which, especially thanks to recent new detections (e.g. Fynbo et al. 
2011; Peroux et al. 2011, 2012; Noterdaeme et al. 2012b; Fynbo 
et al. 2013; Bouche et al. 2013; Jorgenson & Wolfe 2014), has 
yielded 11 DLA galaxies currently known at ^ > 1.9. In Ta- 
ble 2, we present a summary of the properties of these confirmed 
DLA galaxies. This compilation is the result of several heteroge- 
neous searches, which include both serendipitous discoveries, and 
results from programmes that have pre-selected targets according 
to absorption properties, such as high metallicity. For this reason, 
it is quite difficult to reconstruct the selection function of these 
searches, and it is even harder, if not impossible, to obtain a full 
census of the DLAs for which a search has been attempted with- 
out success. Some notable exceptions are the recent VLT/SINFONI 
survey by Peroux et al. (201 1) and Peroux et al. (2012), or the Gem- 
ini/NIFS programme by Wang et al. (MNRAS submitted), as these 
authors report both detections and non-detections. 

A comparison between the properties of the galaxies listed in 
Table 2 and the results of more homogeneous searches, such as our 
study, becomes therefore critical to understand the extent to which 
the confirmed DLA galaxies are representative of the general popu- 
lation of DLA hosts. This comparison is offered in Figure 10, where 
we examine the distributions of impact parameters, metallicity, col- 
umn densities, and SFRs for the known DLA hosts, the galaxies 
detected in our images, and the results of the IFU searches. Based 
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Figure 10. Distribution of impact parameters as a function of SFRs for the 
detected galaxies in proximity of the DLA position (grey crosses), com- 
pared to the distribution of the confirmed DLA galaxies from Table 2 (blue 
squares). Empty squares indicate galaxies without published SFRs, which 
are shown at the arbitrary value of IM© yr“ ^ . The brown triangles mark the 
upper limits on the SFRs for the 9 DLAs without counterparts in the sample 
of Peroux et al. (2012), which we plot at the maximum distance probed by 
their observations. The red solid and dotted lines represent the locus of the 
minimum impact parameters and the corresponding 25 — 75 percentiles (as 
in Figure 7). The black line and upper limits show the maximum impact pa- 
rameter as a function of SFR for a compact LEG that would satisfy the SFR 
surface density inferred for the high-cool population. The green dashed re- 
gion shows instead the 25 — 75 percentiles of the SFR distribution for the 
compact LBGs included in the stack by Rafelski et al. (201 1). 


X 


= 

z 



0.0 

-0.5 

- 1.0 

-1.5 

- 2.0 

-2.5 

22.0 

21.5 
21.0 

20.5 
20.0 



bdia (kpc) 


|X/H] 


Figure 11. Left panels. The metallicity (top) and column density (bottom) 
distributions of known DLA galaxies (blue squares) compared to the distri- 
bution of the sources detected in our survey (grey crosses) for tp ^ 1 M© 
yr“^ and ^ 25%- tiight panel. Same as the central panels, but 

comparing the SFRs in emission with the metallicity in absorption. The lim- 
its from Peroux et al. (2012) and Wang et al. (submitted) are also included 
as downward triangles (filled and open, respectively). SFRs are shown as 
published, and thus these measurements contain additional scatter due to 
inhomogeneous tracers and calibrations. 


on the published star formation indicators for 8/11 DLA galaxies 
(Table 2), it appears that the population traced by the known objects 
is generally above the completeness limit of our survey {tp = 2 Mq 
yr“^). In the following, we will assume that the three DLAs with- 
out published SFRs (Q0139-0824, Q0338-0005, and Q0953+47) 
lie within the sensitivity of our images. We note, however, that the 


compiled SFRs in Table 2 have a substantial uncertainty due to ob- 
servational errors (e.g. slit losses in some cases). They are also quite 
heterogeneous, due to the use of different SFR indicators and cali- 
bration schemes in the literature. 

With these caveats in mind. Figure 10 shows that the known 
DLA galaxies lie all at impact parameters below the typical mini- 
mum impact parameter for a random background population, with 
^ 8/11 systems having impact parameters below bmin, 25 %^- Con- 
versely, our survey indicates that ^ 6/32 of the DLA galaxies 
are found below bmin, 25 % at comparable SFRs. In our survey no 
more than ~ 15% of the DLAs have hosts with > IMq yr“^ at 
low impact parameters. Conversely, previous searches would sug- 
gest instead that these SFRs and impact parameters should arise 
in ~ 70% of the DLAs. This comparison is affected by the small 
sample size, but it highlights how the emission properties of the 
confirmed hosts cannot be generalized to the entire DLA popula- 
tion without considering the various selection effects for the sample 
listed in Table 2. 

Similar considerations can be made from the VLT/SINFONl 
survey by Peroux et al. (2012). This programme provides IFU 
data and hence a complete spectroscopic search up to distances of 
~ 40 kpc in 10 z ^ 2 DLAs with a wide range of absorption 
properties, similar to those shared by our targets. A direct com- 
parison between the SINFONI survey and known DLA galaxies is 
slightly complicated by the wide range of sensitivities reached by 
this programme, with upper limits between i/ ~ 1 — IOMq yr“^. 
Nevertheless, the low detection rate (a single detection in a sam- 
ple of 10 DLAs) reinforces our conclusions (see also Bouche et al. 
2012). Similarly, Wang et al. (submitted) did not find DLA hosts 
in their Gemini/NIFS IFU data down to tp ^ 2.2Mq yr“^ for two 
systems, and tp ^ IIMq yr“^ for a third. 

As discussed, it is far from trivial to quantify the relative im- 
portance of the different selection effects, and here we bring to the 
attention of the reader three of the most obvious ones. First, the 
sample listed in Table 2 suffers from a luminosity bias, which is 
typical of any imaging or spectroscopic search. In this case, the 
confirmed DLA galaxies would represent only the “tip of the ice- 
berg” at the bright end of the luminosity function. The luminosity 
bias alone could explain the excess of sources at small impact pa- 
rameters reported in Figure 10. However, it is worth considering a 
second more subtle selection effect which may arise from the obser- 
vational technique adopted in the discovery of several of the known 
DLA hosts. A non-negligible number of these DLA galaxies have 
been discovered by triangulating the signal detected in two or three 
slits centered at the quasar position (see e.g. figure 1 in Fynbo et al. 
2010; Mpller et al. 2004). In this layout, the likelihood of detect- 
ing galaxies at small impact parameters is higher, purely because 
this is where the area covered by the slits is maximal. While the 
majority of the DLAs would be detected given the model proposed 
by Fynbo et al. (2008), a fraction of DLA hosts may still be unde- 
tected purely because they lie at larger impact parameters. This ef- 
fect would skew the impact parameter distributions of the detected 
DLAs towards small baia, similarly to the observations in Figure 
10 . 

Furthermore, as highlighted in Figure 1 1 , the absorption prop- 
erties of the DLAs with known hosts are skewed towards high 
metallicity (right panel), and slightly towards higher column den- 


® The inequality arises from the fact that the DLA galaxy Q0139— 0824 
at hdla — 13.0 ± 0.4 kpc may intersect the region feuia > 25% for 

v> ~ 1.5 - 7M0 yr-L 
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sities (left bottom panel). The excess of systems at high metallicity 
is a known consequence of the well-documented targeting strat- 
egy. Observations of metal-rich DLAs are indeed favored because 
of a possible mass-metallicity relation (Mpller et al. 2004; Fynbo 
et al. 2008, 2010; Neeleman et al. 2013) that would increase the 
odds of detecting the hosts of the more metal-rich systems. From 
this discussion, we conclude that; i) care should be exerted when 
interpreting the properties of DLA host galaxies solely based on 
compilations like the one in Table 2; ii) models aiming to repro- 
duce DLA galaxies should not only reproduce the scaling relations 
of the known DLA hosts, but also the number of non-detections in 
homogeneous searches (such as the present one). 

Finally, in Figure 11, we also compare the few galaxies de- 
tected in our survey at bjia ^ hmin, 25 % against some previously re- 
ported correlations between absorption and emission properties in 
DLA hosts. Correlations between impact parameters and column 
densities or metallicity have been already reported for the known 
DLAs (e.g. Krogager et al. 2012), as is also shown in the panels of 
Figure 1 1 . Our candidates appear instead to be more scattered, but 
these deviations are not particularly significant in the absence of 
spectroscopic confirmation. When examining how SFRs vary with 
metallicity (right panel of Figure 11), no clear trend is evident in 
either our sample of candidate hosts or in the sample of confirmed 
DLA galaxies. To first order, if a mass-metallicity relation is in 
place for DLAs (Mpller et al. 2004; Fynbo et al. 2008; Neeleman 
et al. 2013; Mpller et al. 2013), one would expect a similar correla- 
tion between SFRs and metallicity, given the known correlation be- 
tween SFRs and masses in galaxies. Flowever, due to the small sam- 
ple size and the intrinsic errors in the SFRs mentioned above, this 
lack of correlation is inconclusive. This is especially true for our 
candidates without redshift confirmations. Furthermore, the trend 
of metallicity with impact parameter complicates the relation be- 
tween metallicity in absorption and SFR integrated over the entire 
host galaxy (Fynbo et al. 2008). It would nevertheless be interesting 
to augment the size of uniformly-selected DLA hosts to investigate 
further how a mass-metallicity relation for DLAs is reflected into a 
metallicity-SFR correlation. 

3.4 Implications for cooling rates and emission from LBG 
outskirts 

Having characterized the distribution of galaxies in the suiTound- 
ings of DLAs, we are now able to comment on the last possible 
scenario for the heating mechanism of the high-cool DLAs, which 
we discussed in Section 2.4. By means of the limits on the in-situ 
SFRs, we highlighted a potential discrepancy between the SFRs 
that would provide the heating rate necessary to power high-cool 
DLAs and the limits inferred by our direct observations for both 
extended low-surface brightness sources and compact (r < 0.25”) 
LBGs embedded in the DLA gas at very small impact parameters. 
As noted, this tension could be resolved if the heat is provided by 
LBGs at larger impact parameters, an hypothesis we can now ad- 
dress explicitly. 

Due to the r~^ dependence of the UV flux, LBGs at larger dis- 
tances have to be much more luminous to provide the same photon 
flux at the DLA position. In Figure 10, we compare the locus of im- 
pact parameters as a function of SFRs for all the galaxies detected 
near the targeted DLAs with the maximum impact parameters at 
which we can find an LBG that satisfies the constraints imposed by 
the DLA cooling rates (cf. Figure 6). This limit is computed as fol- 
lows. Given the size of LBGs ribg, the size of the associated DLA 
rdia (and hence of the maximum impact parameter) is constrained 


by the area covered by DLAs relative to the area covered by LBGs, 
C'dia.ibg. Wolfe & Chen (2006) and Wolfe et al. (2008) estimate 
Cdia.ibg ~ 330 between z ~ 2.5 — 3.5. We can therefore com- 
pute the SFR required to satisfy the SFR surface density inferred 
for the high-cool population as = EsnTrrdia = EsfrvrrfbgChc, 
where the covering fraction of the high-cool population is Che ~ 
l/2Cdia,ibg, as discussed in Wolfe et al. (2008). 

Figure 10 shows the resulting limits on the impact parameters 
for LBGs, where the minimum at ~ 2 M 0 yr“^ is imposed by 
the constraint ribg ^ 0.5 kpc. In this calculation we further as- 
sume Esfr ~ 10“^ ® M 0 yr”’^ kpc“^ for the high-cool DLAs from 
the median of the distribution shown in Figure 6. We note that the 
presence of a dust-obscured population, unexplored in our study, 
would not contribute towards the budget of UV photons needed to 
satisfy the model prediction for the high-cool DLAs, as these pho- 
tons will be screened by the same dust that prevents the detection 
of the galaxies in our imaging. 

Given our sample of 32 DLAs and assuming a roughly equal 
split between low-cool and high-cool DLAs, we would expect to 
have detected approximately 8 galaxies below the impact parameter 
limits imposed by the calculation above. Our data show that ^ 6 
candidate DLA hosts (without double-counting the two detections 
in 5:G5 field) are found below this line. If we allow for possible 
leakage in the 10;G11 field, only ^ 5 galaxies are found out of 
the ~ 8 galaxies expected. Currently, these numbers are in formal 
agreement, but the modest number of star-forming galaxies below 
~ 30 kpc may become uncomfortable for the model of Wolfe et 
al. (2008) if some of these candidates were found to not be at the 
DLA redshift. However, any tension should be considered modest, 
if present at all, given the many uncertainties at play. Firstly, this 
comparison hinges on a small number of galaxies; secondly, the 
fraction of high-cool DLAs relative to low-cool DLAs is unknown 
in our sample and may differ from what has been reported by Wolfe 
et al. (2008); thirdly, previous estimates of the relative covering 
fraction of DLAs and LBGs rely on shallower luminosity functions 
than those now available (e.g. Alavi et al. 2014). It will be very 
interesting to use larger samples and new direct measurements to 
investigate the assumptions at the base of the model by Wolfe et al. 
(2003b,a). 

Finally, with a census of the LBGs in the surroundings of these 
DLAs, we can compare our findings with the model proposed by 
Rafelski et al. (2011), according to which DLAs arise in the out- 
skirts of compact LBGs. In Figure 10, we show the 0.25 — 0.75 
percentiles of the SFRs for the galaxies included in the stack of 
Rafelski et al. (201 1), which lie in the range ip = 0.8 — 2 M 0 yr“^. 
Unfortunately, our survey starts being incomplete below ip = 2 M 0 
yr“^, preventing a statistical comparison between our detections 
and the presence of compact LBGs. Again, deeper imaging surveys 
will provide valuable constraints for this model. 


4 WHAT ARE DLAs? 

A definite answer to the question of what DLAs are has been hard 
to obtain, despite much progress in characterizing the emission 
and absorption properties of these absorbers. Besides the interest 
in unveiling the nature of these systems, a deeper understanding 
of DLAs has far reaching implications. DLAs represent the major 
repository of neutral gas in the « ~ 2 — 3 Universe, and studies of 
the SFRs in representative samples of these absorbers, like the one 
presented here, offer new insight into how the most significant HI 
overdensities relate to star formation. 
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The primary results emerging from our analysis are: 1) new 
limits on the in-situ SFR of DLAs, to i/t < 0.09 Mq yr“^ for com- 
pact regions or t/> < 0.27 Mq yr“^ for more extended sources; 
2) a deep search for galaxy counterparts in the surroundings of 
DLAs, through which we have detected an excess of sources only 
at small impact parameters (< 30 kpc), and for < 4/32 DLAs 
to a completeness limit of = 2 Mq yr“^. Thus, our survey 
places new quantitative limits on the allowed parameter space for 
DLA galaxies, ruling out the parameter space tp ^ 2 M 0 yr“^ and 
hdia < &min, 25 % kpc for the general DLA galaxy population. 

In this section, we aim to discuss how these new limits fit with 
the empirical knowledge accumulated over the years on the origin 
of DLAs. In particular, we will focus on three scenarios: i) DLAs as 
systems where star formation occurs only throughout the absorbing 
gas; ii) DLAs as moderately star-forming galaxies which are how- 
ever sufficiently dusty to become invisible at UV wavelengths; iii) 
DLAs as UV-bright star-forming galaxies, or LBGs. At this stage, 
the discussion will be qualitative and, occasionally, even specu- 
lative. Future theoretical work is encouraged to incorporate more 
quantitatively the new empirical constraints we have provided in 
the context of the many models of DLAs that can be found in the 
literature (e.g. Fynbo et al. 1999; Nagamine et al. 2007; Fynbo et 
al. 2008; Pontzen et al. 2008; Tescari et al. 2009; Hong et al. 2010; 
Cen 2012; Rahmati & Schaye 2014; Bird et al. 2014; Barnes & 
Haehnelt 2014; Danovich et al. 2014). 

Furthermore, in the following, we will not explore the even 
harder question of the moiphology of DLAs, which can take the 
form of disk-like structures, clumps, tidal debris, inflows, outflows, 
or even quite extended structures (Giovanelli & Haynes 1989). We 
will simply focus on the more general relationship between DLA 
gas and star formation. It is also implicit that a rigorous classifica- 
tion within the three classes mentioned above is not mandatory, or 
may even be incorrect if DLAs arise from a diverse population. 

4.1 DLAs as aggregates of gas and stars 

In this first class, we envision DLAs as aggregates of Hi gas giv- 
ing rise to the absorption signature, with star formation occur- 
ring throughout the same gas. In this scenario, the DLA gas and 
the DLA galaxy are coincident in the same system, which can be 
probed either in absorption (the DLA gas) or in emission (the DLA 
galaxy). The plausibility of this scenario has already been investi- 
gated in the literature, most notably by Wolfe & Chen (2006), or 
more recently by means of composite spectra of DLAs (Rahmani 
et al. 2010; Noterdaeme et al. 2014; Cai et al. 2014). In agreement 
with these previous studies, our direct measurements of the in-situ 
SFRs via FUV emission rule out that typical DLAs are star-forming 
pockets of gas that are buried under the glare of the background 
quasars, up to SFR surface density limits of Ssfr = 10“^ ®'^ Mq 
yr“^ kpc“^ and Sgfr = Mq yr“^ kpc“^, for sizes of 

2 kpc and 12 kpc, respectively (Figure 4). 

During our analysis, we emphasized that the low in-situ SFRs 
inferred from our observations are not inconsistent with other ob- 
served properties of DLAs, essentially allowing for a DLA popula- 
tion with typically very little star formation. As shown in Figure 5, 
even modest star formation events can still satisfy the local metal 
enrichment. Additionally, provided that most of the DLA gas at 
higher redshift merges onto other structures that are forming stars 
at higher rates (Prochaska & Wolfe 2009), we can still preserve the 
general argument that DLAs contain enough neutral hydrogen to 
form about a third of the stars produced in the universe by 2 ; ~ 0, 
even without imposing that these stars are formed in situ. 


Therefore, our observations do not exclude the possibility that 
star formation occurs throughout the DLA gas at rates below our 
inferred limits, leaving open a scenario in which stars are formed 
fully within the DLA gas for the origin of most DLAs. However, 
as argued in Section 2.2, this scenario is more difficult to corrob- 
orate from a theoretical point of view: typical DLAs with column 
densities Nai = — 10^^ cm“^ may not provide a favor- 

able environment for star formation. Direct observations are needed 
to test this hypothesis, but we note that, in line with our argu- 
ment, there is evidence for in-situ star formation only in DLAs with 
Nm > 10^^ cm“^ (Rafelski et al. 2011; Noterdaeme et al. 2014), 
and not at lower column densities. Furthermore, the potential dis- 
crepancy between the observed limits and the inferred SFR surface 
densities from the [C ll*] cooling rates (Figure 6) may be resolved 
by associating “high-cool” DLAs to nearby star-forming galaxies 
(Wolfe et al. 2008), a fact that would make a picture in which DLAs 
form stars fully within the absorbing gas insufficient to account for 
all the known observables. 

And while recent studies have uncovered a population of “dark 
galaxies” (Rauch et al. 2008; Cantalupo et al. 2012), despite the evi- 
dence that DLAs lack appreciable star formation in situ, we believe 
that it would be incorrect to extrapolate our findings to conclude 
that DLAs simply arise from systems not related to star formation. 
In fact, the general DLA population is not an aggregate of pristine 
gas, nor does it share the same metal content of the IGM (e.g. Rafel- 
ski et al. 2012). Furthermore, the median velocity widths of DLAs 
(Prochaska & Wolfe 1997) suggest that DLA gas resides in signif- 
icant potential wells, in which it is plausible to expect star-forming 
regions. 

We conclude that a scenario in which DLAs are aggregates of 
gas and stars is not completely ruled out by observations, although 
it is not the most favorable hypothesis, given the aforementioned 
issues. Thus, if not fully embedded with DLA gas, star formation 
has to occur at least in the proximity of DLAs, a scenario we will 
explore in the next two sections. 

4.2 DLAs as dusty star-forming galaxies 

In the two scenarios we discuss in this and in the following section, 
we consider DLAs that are associated with some star formation, 
which is however not fully embedded in the absorbing gas, as it 
was instead the case for the scenario in Section 4.1. The readers 
who are familiar with the work by Wolfe & Chen (2006) and Wolfe 
et al. (2008) will recognize this as the “bulge hypothesis”. 

A substantial fraction of the total cosmic star formation arises 
from obscured regions (e.g. Blain et al. 1999; Chary & Elbaz 2001; 
Bouwens et al. 2009; Decarli et al. 2014), and given that DLAs 
contain enough neutral hydrogen to account for about a third of the 
stars ever formed, we consider whether most DLAs may originate 
in proximity to dusty star-forming regions. In this category, we do 
not envision associations between DLAs and submillimetre or ultra 
luminous infrared galaxies for a few reasons. Firstly, the space den- 
sity of these extreme populations is insufficient to account for the 
more abundant DLA population. Secondly, the few direct searches 
in emission for molecular gas in DLAs rule out the presence of 
massive reservoirs of cold gas and dust near the absorbing gas (e.g. 
Wiklind & Combes 1994, , Dessauges-Zavadsky et al. in prep.). 

Thus, with the expression “dusty star-forming galaxies”, we 
identify instead a population of galaxies with masses comparable to 
or even lower than the classic LBGs (Mhaio ~ — 10^^ Mq), 

in which star formation occurs within regions that are dusty enough 
to be attenuated and become invisible to surveys that select galaxies 
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solely based on FUV emission. Several studies (e.g. Reddy & Stei- 
del 2004; Reddy et al. 2010) concluded that dust in bright LBGs 
suppresses the FUV emission by a factor of ~ 4 — 5. Extinction up 
to ~ 2 mag, if extrapolated also to fainter systems, could therefore 
play an important role in shaping the statistics of the impact pa- 
rameter as a function of SFR (Figure 7), especially given that our 
results are based on rest-frame wavelengths < 1000 — llOOA. 

Therefore, the presence of a relatively massive but dusty 
galaxy population near DLAs should be considered further, as it 
would for instance ease the potential tension between studies that 
reveal how DLAs originate from regions with high bias (Cooke et 
al. 2006; Font-Ribera et al. 2012) and theoretical works that instead 
struggle to incorporate both this high bias and all the other observed 
DLA properties in a single model (e.g. Bird et al. 2014; Barnes & 
Flaehnelt 2014). 

Flowever, measurements of the in-situ metallicity (e.g. Figure 
5) place DLA gas in regions of modest dust content. This is also 
seen directly in the low reddening measurements along the line 
of sight to quasars with DLAs (Vladilo et al. 2008; Khare et al. 
2012), although we cannot fully rule out with current observations 
the presence of dusty clumps in close proximity to, but not overlap- 
ping with, the locations probed in absorption. Similarly, the paucity 
of cold gas and molecules in high-redshift DLAs (see e.g. Carilli et 
al. 1996; Kanekar et al. 2006; Curran & Webb 2006; Jorgenson et 
al. 2006; Noterdaeme et al. 2008; Ellison et al. 2012; Srianand et 
al. 2012; Kanekar et al. 2014; Jorgenson et al. 2014) points towards 
an environment that is atypical for dusty, gas rich, galaxies. Finally, 
we note that, although the exact relationship between the observed 
cooling rate and the UV photon flux may be subject of future revi- 
sion, the presence of high-cool DLAs disfavor a scenarios in which 
the UV radiation field is completely suppressed by dust. 

In summary, given these contrasting pieces of evidence, the 
possibility that dust obscuration plays a role in shaping the detec- 
tion rates of DLA galaxies remains open. This scenario is currently 
almost completely unexplored for galaxies of modest IR luminos- 
ity, and unbiased surveys at infrared and millimeter wavelengths 
to very faint levels are now needed to address the importance of 
dusty galaxies of modest mass as an important or even dominant 
population of DLA hosts. 

4.3 DLAs as UV-selected star-forming galaxies 

As already noted, the absorption properties of DLAs (e.g. metallic- 
ity and cooling rates), the fact that DLAs contain most of the neutral 
gas reservoir at high redshift, and the results of cross-correlations 
of DLAs and LBGs suggest that DLAs arise in the vicinity of star 
forming galaxies. With the caveat of a possible dusty population 
discussed in Section 4.2, it is therefore natural to associate DLAs 
with LBGs, as done by many authors in the literature (e.g. Fynbo et 
al. 1999; Wolfe et al. 2003b; Mpller et al. 2004; Fynbo et al. 2008; 
Pontzen et al. 2008; Rafelski et al. 2011; Rahmati & Schaye 2014). 
Given that DLAs are selected purely based on gas cross-section and 
not luminosity, and because of the difficulties of identifying DLA 
host galaxies, it seems natural to associate DLAs with the faint end 
of the luminosity function (Fynbo et al. 2008; Rauch et al. 2008; 
Rafelski et al. 201 1). In this scenario, which is also corroborated by 
simulations (Haehnelt et al. 2000; Nagamine et al. 2007; Pontzen 
et al. 2008; Rahmati & Schaye 2014), the few known DLA hosts 
with i/i > 1 — 2 Mq yr“^ arise from the brighter LBGs. 

Our results are consistent with this picture, providing for the 
first time quantitative limits on the bright end of the luminosity 
function of the DLA host galaxies. Our observations rule out the 


possibility that the DLA galaxies are LBGs with ip ^ 2 Mq yr“^ 
(unobscured) which are outshined by the background quasars. Fur- 
thermore, the lack of LBGs with ip ^ 2 Mq yr“^ and baia ^ 
^>min, 25 % kpc (Figure 7) rules out “bright” LBGs with large Hi 
disks or more amorphous structures as those predicted by modem 
simulations (e.g. Danovich et al. 2014) of sizes ~ 5 — 20 kpc as 
the dominant population that gives rise to DLAs. At face value, this 
empirical fact disfavors models in which most of DLAs originate 
from extended rotating disks (Prochaska & Wolfe 1997; Mailer et 
al. 2001). 

Given these considerations, we are therefore forced to con- 
clude that DLAs arise in the surroundings of fainter galaxies, with 
typical SFRs of ip <2 Mq yr“^ (e.g. Rauch et al. 2008; Rafelski et 
al. 2011). If these dwarf LBGs are centrals of their dark matter ha- 
los (i.e. if they are field galaxies), then it would entirely explain the 
lack of bright LBGs in the surroundings of DLAs (e.g. Figure 9). 
However, a scenario in which DLAs are typically isolated dwarf 
LBGs may be difficult to reconcile with the measured DLA bias 
(Cooke et al. 2006; Font-Ribera et al. 2012), which suggests instead 
a comparable clustering amplitude of DLAs and bright LBGs. As 
speculated by Font-Ribera et al. (2012), this tension may be allevi- 
ated if DLAs arise from dwarf galaxies which are instead satellites 
of more massive (brighter) LBGs, but a more quantitative analysis 
is now needed. Incidentally, we note that the idea of an extended 
DLA gas shaicture, in which sub-halos are embedded and which is 
larger than the typical H I disk, is what motivated the simple picture 
of extended disks put forward by early models (Prochaska & Wolfe 
1997; Mailer et al. 2001) to explain the observed kinematics. And 
this picture gains also support from modern cosmological simula- 
tions, in which DLA gas extends inside or in proximity to massive 
dark matter halos, also encompassing smaller sub-halos (Rahmati 
& Schaye 2014). 

To muddle the picture, however, simulations still struggle to 
reproduce correctly some of the most basic properties of DLAs, 
such as the redshift evolution of their number (e.g. Bird et al. 2014). 
Also, the low rate of detections in the IFU observations of Peroux 
et al. (2012) and Peroux et al. (201 1) out to 40 kpc from the DLA 
position may not be easy to reconcile with this idea, as one would 
expect to detect the bright LBGs at the center of the massive parent 
halos at these distances in some instances. Perhaps more critical is 
the lack of significant excess in the DLA-LBG correlation function 
beyond 20 kpc, which is at odds with the idea of an abundant DLA 
population in satellites of massive LBGs out to their virial radii 
(Figure 9). 

We therefore conclude that, while the idea of a link between 
LBGs and DLAs is sound from a qualitative point of view, more 
work is needed to quantitatively incorporate the properties of the 
observed DLA galaxies (e.g. Krogager et al. 2012), the limits of 
the bright end of the luminosity function (Figure 7), the DLA bias, 
kinematics, and the presence of metals in high-ionization states 
(Lehner et al. 2014) in a coherent model. Deep and complete red- 
shift surveys in the surroundings of these DLAs are further needed 
to better constrain the importance of clustering of dwarf galaxies 
around massive LBGs in shaping the properties of the DLA host 
galaxies. These searches should however account for the subtle bi- 
ases induced by the stochastic nature of star formation at these low 
masses, which shapes the detection rates of faint galaxies (da Silva 
et al. 2014; Dominguez et al. 2014). 
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5 SUMMARY AND CONCLUSIONS 

In this paper, we have presented results from an imaging survey of 
32 quasar fields hosting DLAs at 2 : ~ 1.9 — 3.8, with a median 
redshift z ~ 2.7. By leveraging the double-DLA technique that 
allows us to image the FUV continuum of DLA host galaxies at 
all impact parameters, including the innermost regions traditionally 
outshined by the background quasars, we have directly studied the 
in-situ emission properties of the DLA gas and the distribution of 
the surrounding galaxies. We have also systematically compared 
our findings to the properties of DLAs inferred from absorption 
spectroscopy. Our main results are summarized below. 


- By positioning apertures of 2 kpc and 12 kpc centered on the 
DLA in each field, we have constrained the median in-situ SFRs of 
the DLA gas to be t/i < 0.38 Mq yr“^ and ip < 0.65 M© yr“^ in 
individual systems. With a stacking analysis, we have derived even 
tighter constraints on the median SFRs of DLAs, ip < 0.09 Mq 
yr“^ and ip < 0.27 Mq yr“^ in the 2 kpc and 12 kpc aper- 
tures. Converted into SFR surface densities, these limits become 
logSsfr < —1.54 M 0 yr“^ kpc“^ and logSsfr < —2.62 Mq 
yr“^ kpc“^, ruling out the existence of appreciable star formation 
embedded in the DLA gas. These limits reinforce the statistical 
conclusion of Wolfe & Chen (2006) and are in line with expec- 
tations borne from the current theoretical understanding of the star 
formation law. 

- By comparing the limits on the in-situ star formation to the 
observed metal content of DLAs in absorption, we have found that 
DLA gas can be enriched to the observed levels in > 0.05 — 10 
Myr for compact DLAs and in > 1 — 300 Myr for extended DLAs. 
Given the available time for the enrichment to occur (~ 2 Gyr), 
even lower rates of star formation (30 — 500 times lower than the 
measured upper limits) would suffice to enrich DLAs to the ob- 
served values, without the need for external sources of metals. It 
thus appears unlikely that star formation in or near DLAs underpro- 
duces the metals locked in DLAs, in line with the fact that DLAs 
contain only a modest fraction (~ 1%) of all the metals produced 
in high redshift LBGs. 

- When compared to the inferred SFR surface densities required 
to satisfy the [C ll*] cooling rates of the “high-cool” DLA popula- 
tion, we find a tentative discrepancy between our measured lim- 
its and the star formation required to heat the DLA gas. Further- 
more, the number of star-forming galaxies in the surroundings of 
the DLAs with our deep imaging is slightly lower than the num- 
ber expected to satisfy the heating rates predicted by models. This 
tension is however only very marginal and should be investigated 
further in larger samples. 

- By studying the distribution of impact parameters as a func- 
tion of the star formation rates for all the galaxies detected in prox- 
imity to the DLA position, we have identified 6 possible fields 
with an excess of sources compared to expectations from a random 
galaxy background. Additional follow-up spectroscopy is needed 
to confirm the redshift of these candidate DLA galaxies. 

- At the completeness limit of our survey (ip = 2Mq yr“ ^ ), we 
find that ^ 13% of the DLAs lie at impact parameters feaia/kpc ^ 
{iP/Mq yr“^)° ® -|- 6, in contrast with samples of known DLA 
hosts which appear to suffer from selection effects in luminos- 
ity and/or impact parameters. Furthermore, we do not detect ex- 
cess in the projected DLA-galaxy correlation function between 
20— 100 kpc, ruling out an abundant population (2> 50%) of bright 
DLA host galaxies at these separations. 


Our study adds new direct constraints to the debate of what 
DLA galaxies are, albeit not solving this puzzle. Our new find- 
ings, combined with findings from other studies, suggest that DLAs 
cannot originate from aggregates of gas which are forming stars 
throughout to levels of ~ 0.1 — 0.6 Mq yr“^. Through indirect ob- 
servations and theoretical arguments, we argued that a more plau- 
sible scenario, which has been already advocated for by many au- 
thors, puts DLA gas in proximity to star-forming galaxies, possibly 
UV-selected LBGs, although the presence of a bias induced by dust 
needs to be investigated further. While it appears evident that DLAs 
are not directly associated to bright LBGs and that dwarf galaxies 
are the more natural hosts, one could envision associations with ei- 
ther dwarf isolated galaxies, or dwarf galaxies which are clustered 
with more massive LBGs. The latter case may be required to re- 
produce the observed bias of DLAs, but it appears to be at odds 
with the paucity of bright LBGs found in previous IFU searches 
or as suggested by our measurement of the DLA-galaxy projected 
correlation function. Further modeling is now required to quantita- 
tively explore how to incorporate our new constraints with earlier 
observations and models for the origin of DLAs. 

This study further opens new prospects to finally unveil the na- 
ture of DLAs. Our survey has placed strong limits on the presence 
of star formation in galaxies close to the quasar sightline. An im- 
portant open question that future redshift surveys should address is 
the distribution of star-forming galaxies at large impact parameters 
to the sightline. Large integral field spectrographs, such as MUSE 
at VLT and KWCI at Keck, will soon provide this piece of infor- 
mation. Similarly, the deployment of the full ALMA array should 
enable the first glimpse into a potential dust obscured galaxy pop- 
ulation invisible to UV studies, which may be associated with (a 
fraction of) the DLA population, but which would have been com- 
pletely missed by our study. 

Furthermore, we highlighted on several occasions the poten- 
tial of a deeper imaging survey which adopts the double-DLA tech- 
nique at the basis of our study. A similar survey, but with an im- 
provement by a factor of 10 in depth, would yield stringent con- 
straints on the in-situ properties of DLAs, critical to investigate the 
SF law for low metallicity gas and the mechanisms that are respon- 
sible for the enrichment and heating of the largest neutral hydrogen 
reservoir in the high-redshift Universe. Furthermore, these obser- 
vations would provide unprecedented information on the clustering 
of faint dwarf galaxies around DLAs, as as well as a direct discrim- 
inant between models for DLAs. Thus, carefully designed surveys 
at current facilities and, in the future, observations at the 30 m tele- 
scopes have the potential of finally answering some of the most 
challenging questions that have characterized DLA studies in the 
past thirty years. 
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